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Introduction 

In  a  previous  paper  [l]1  an  attempt  was  made  to  contribute 
toward  a  comprehensive  and  —  hopefully  —  rigorous  analytical 
theory  of  linear  viscoelastic  behavior,  comparable  in  scope  and 
character  to  the  systematic  body  of  general  propositions  that  has 
long  been  available  in  the  classical  theory  of  elastic  solids. 

The  material  contained  in  [l],  which  is  confined  to  an  isothermal 
quasi-static  treatment  of  the  subject,  includes:  results  concern¬ 
ing  the  structure  of  the  relevant  constitutive  relations  and  the 
connection  between  various  alternative  versions  of  the  stress- 
strain  law;  conclusions  regarding  the  nature  of  the  time  and 
position  dependence  of  solutions  to  the  governing  field  equations; 
integral,  reciprocal,  and  uniqueness  theorems;  and  theorems  per¬ 
taining  to  the  integration  of  the  field  equations  in  terms  of 
complete  systems  of  stress  functions.  Later  on,  Gurtin  [2] 
obtained  generalizations  to  viscoelasticity  theory  of  the  known 
variational  principles  of  elastostatics .  Also,  certain  results 
appearing  in  [l]  have  since  been  extended  to  non-isothermal  condi¬ 
tions  in  [33,  [4]. 

In  the  present  paper  we  continue  the  project  begun  in  [ 1 ] 
and  turn  to  the  study  of  Green's  functions  in  the  isothermal 
quasi-static  theory  of  a  homogeneous  and  isotropic  viscoelastic 
medium  which  obeys  the  general  linear  relaxation  integral  law. 

1  Numbers  in  brackets  refer  to  the  list  of  publications  at  the 
end  of  this  paper. 
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Our  main  objective  in  this  connection  is  twofold:  first,  we  aim 
at  integral  representations  for  the  solution  of  the  standard 
boundary- value  problems,  analogous  to  those  originated  in 
elasticity  theory  by  Betti,  Volterra,  Lauricella,  and  Somigliana; 
second,  we  seek  to  apply  such  integral  representations  to  a  proof 
of  Saint-Venant 1 s  principle  in  a  formulation  that  is  the  counter¬ 
part  in  viscoelasticity  of  the  elastostatic  principle  suggested 
by  von  Mises  [6]  and  established  in  [7]. 

The  chief  mathematical  tool  in  this  investigation  is 
once  again  supplied  by  the  algebra  and  calculus  of  Stieltjes 
convolutions  developed  in  [l],  Further,  some  of  the  results  on 
linear  viscoelasticity  appearing  in  [l]  are  essential  prerequisites 
to  our  current  purpose.  For  this  reason,  and  in  order  to  render 
the  present  paper  sensibly  self-contained,  we  collect  in  Section  1 
various  definitions  and  theorems,  mostly  adapted  from  [l],  which 
are  needed  repeatedly  in  the  subsequent  analysis. 

In  Section  2  we  deal  first  with  the  problem  of  a  con¬ 
centrated  load  acting  at  a  point  of  a  viscoelastic  solid  that 
occupies  the  entire  space.  The  solution  to  this  singular  problem, 
which  is  a  generalization  of  Kelvin's  problem  in  elastostatics, 
is  defined  and  deduced  explicitly  by  means  of  a  limit  process 
applied  to  the  solution  of  a  sequence  of  regular  problems  governed 
by  distributed  body  forces .  The  basic  singular  solution  of  the 


2 


See  Love  1 5 3 »  Art.  169,  for  detailed  references. 
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field  equations  thus  established  is  subsequently  used  to  generate 
the  higher-order  singularities  appropriate  to  force  doublets  with 
and  without  moment.  Certain  relevant  properties  of  the  singular 
solutions  arrived  at  in  this  section  are  studied  in  detail . 

The  results  obtained  in  Section  2  are  applied  in  Section  3 
to  the  contraction  of  Green's  functions  and  the  derivation  of 
integral  representations  for  the  solution  to  the  fundamental 
boundary- value  problems  in  linear  viscoelasticity  theory.  Both 
Section  2  and  Section  3  are  influenced  by  a  partly  parallel 
treatment  in  [8]  of  the  analogous  topics  in  the  classical  equi¬ 
librium  theory  of  elastic  solids. 

The  integral  representations  obtained  in  Section  3  are, 
in  turn,  employed  in  Section  4  to  prove  a  Saint-Venant  principle 
appropriate  to  viscoelastic  solids  within  the  theoretical  frame¬ 
work  underlying  this  paper. 
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1.  Notation .  Preliminaries . 

In  this  section  we  cite  —  occasionally  in  a  modified  or 
extended  form  adapted  to  our  present  needs  —  certain  definitions 
and  results  appearing  in  [1]  that  will  be  required  repeatedly 
later  on. 

Throughout  this  paper  the  letter  E  designates  a  three- 
dimensional  Euclidean  space.  Ihe  symbol  ft,  in  the  absence  of 
any  qualifying  restrictions,  will  denote  an  arbitrary  region  in 
E  that  may  be  either  open  or  closed.  Further,  we  employ  the 
standard  notation  for  closure;  thus,  if  ft  is  open,  stands  for 
the  closure  of  ft.  On  the  other  hand,  the  letter  R  is  consist¬ 
ently  reserved  for  a  regular  region  in  E,  by  which  we  mean  an 
open  (not  necessarily  bounded  or  simply  connected)  region,  the 
boundary  B  of  which  consists  of  a  finite  number  of  non-intersect¬ 
ing  closed  regular  surfaces,  the  latter  term  being  used  in  the 
sense  of  Kellogg  l9)(p.H2).  Since  B  is  bounded  even  if  R  is 
not,  an  unbounded  R  is  of  necessity  an  exterior  region,  i.e.  a 
region  containing  all  sufficiently  distant  points.  Note  also 
that  B  may  have  corners  and  edges.  A  point  P  of  B  will  be 
referred  to  as  a  regular  boundary  point  if  B  possesses  a  tangent 
plane  at  P;  by  a  regular  subset  of  B  we  shall  mean  one  that 
consists  exclusively  of  regular  boundary  points. 

In  the  present  investigation  we  have  frequent  occasion 
to  deal  with  real -valued  functions  of  position  and  time,  whose 
domain  of  definition  is  the  cartesian  product  of  a  region  of 
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apace  6  and  an  (open,  closed,  or  half-open)  Interval  of  time, 
for  which  we  use  the  symbol  3.  In  this  physical  context  we  shall 
denote  by  x,  with  the  rectangular  coordinates  (x^x^x^),  the 
position  vector  of  points  in  R,  call  t  the  time,  and  write  8*3 
for  the  appropriate  domain  of  definition.  Further,  if  f  is  a 
function  defined  on  ftxO,  we  write  f (x,t)  for  the  value  of  f  at 
(x,t)  and  use  f(*,t)  to  designate  the  subsidiary  mapping  of 
position  that  results  from  holding  t  fixed  in  3*  The  analogous 
interpretation  applies  to  f(x, •)•  Finally,  as  far  as  the  partial 
space  and  time  differentiation  of  f  is  concerned,  we  adopt  the 
notation 


.(n) 

,1J 


d’^ffx.t) 


.(x,t) 


m  indices 


dxjdxj . dx^dt11 


d^ffx.t) 


(m,n  =  0,1,2,  •••  •)  , 


(1.1) 


8  35yKJ....:ax^t  (»-  0.1,2,  ••••)  ,  J 

m  indices 

with  the  understanding  that  all  subscripts,  unless  otherwise 
specified,  henceforth  have  the  range  of  the  integers  (1,2,3). 

Turning  to  functions  of  the  time  alone,  we  define  the 
Heaviside  unit  step  function  by 


(a)  h(t)  =  0  for  t£(-oo,0), 

(b)  h(t)  *  1  for  t£[0,oo)  , 


J 


(1.2) 


and  introduce  the  following  convenient  generalization  of  the 
Heaviside  function. 
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Definition  1.1  (Functions  in  Heaviside  Class  H^).  We  say  that 
f  €  H1*  if  f(t)  is,  defined  for  all  t  6  (-00, oo  )  and 


(a)  f  =  0  on  (-oo,0), 

(b)  f  S  CN([0,oo  )). 

Further,  if  f  €  H*1,  we  write 

f(n)(0)  e  f^n^(Of)  (n  -  0,1,2, ...N).  (1.4) 


A  useful  extension  of  this  definition  to  functions  of 
position  and  time  is  furnished  by 

Definition  1.2  (Functions  in  Class  CM,N  or  Class  HM,N) . 

(A)  We  say  that  f  £  CM,N(ftx3)  if  f(x,t)  1b  defined  for 
all  (x,t)£ftxQ  and  the  functions 


r(n) 

,ij . k  (m  =  0,1,2, ---M;  n  -  0,1,2,  •  • -N)  (1.5) 

v - - - > 

m  indices 

exist  and  are  continuous  on  £*3. 

(B)  We  say  that  f 6  HM,N(ft)  if  f (x,t)  is  defined  for 
all  (x,t)  €  Wx(-co  ,00 )  and 


(a)  f  =»  0  on  (2x(-oo,0), 

(b)  f  £CM'N(«x[o,co)). 
Further,  if  f£HM'N(k),  we  write 


.(n) 

,ij 


.(*,0) 


.(n) 


m  indices 


m  indices 


.(*,0+)  (m 


n 


0,1,2, ---M; 
0,1,2, ••*N). 


(1.7) 


3  If  X  is  a  set  then  C^X)  stands  for  the  set  of  all  functions 
that  are  N  times  continuously  differentiable  on  X. 
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Next  we  collect  some  basic  results  from  the  theory  of 
Stieltjes  convolutions.  To  this  end  suppose  that  9  and  i|>  are 
functions  of  position  and  time  defined  on  ®x[o,oo  )  and  fc*(-oo  ,00 ), 
respectively,  and  assume  that  the  Riemann-Stieltjes  integral 

i7(x,t)  =  J  9(x,t-x)d4i(x,x)  (1.8) 

T=-00 

exists  for  all  (x,t)  €  ft  x(-oo  ,00  ) .  Then  the  function  so 
defined  on  ftx(-oo,co)  is  said  to  be  the  Stieltjes  convolution  of 
9  and  we  also  write 

•0's  9#d4>  ,  Otx,t)  »  [9-»ci»|»  J(x,  t) .  (1.9) 

In  view  of  Theorem  1.2  and  Theorem  1.6  of  [1],  we  have 
Theorem  1.1  (Properties  of  Stieltjes  convolutions).  Let 
9€HM'1(G)  and  i|>,©€  H*1'1^) .  Then: 


(a) 

9*a|>€HK»1(^)  with  K  -  min(M,N); 

(b) 

9#dij>  =  \J)*d9  j 

(c) 

9*d('t,*d0)  =  (9#dv|>)«d©  =  9»d9*dQ ; 

(d) 

9*d(i|H-©)  *  9*d>|>  +  9*d©; 

(e) 

9#dtJ>  «=  0  — »  9  =  0  OT  = 

=  0; 

(f) 

9*dh  =  9; 

t 

(g) 

[9*dij>3(x,t)  =  ij>(x,0)9(x,t) 

+  J  9(x,t-x)ij»(T)dT 

for  all  (x,t)€  i.  x(o,oo  ) ; 

0 

(h) 

(i|>*d©)  .  =  tji  ,#d©  +  9*d®  . 

,1  ,1  j  1 

if  N  >  1. 

Theorem  1.3  and  Theorem  1.4  of  [1]  may  be  combined  into 
Theorem  1.2  (Stieltjes  inverse).  Let  9CH2  and  9(0)/0.  Then 
there  exists  a  unique  f miction  \|>  &  9"^,  which  we  call  the 
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Stlelt,1e8  Inverse  of  9,  such  that  \J»€  H1  and 

9*d^  =  h  on  (-00,00).  (1.10) 

Finally,  we  state  and  prove  a  result  on  Stieltjes  convolu¬ 
tions  that  is  closely  related  to  Theorem  1.5  in  [l]. 

Theorem  1.3  (Sequences  of  Stieltjes  convolutions).  Let  ft  be 
bounded  and  let  [911}  be  a  sequence  of  functions  such  that 

(a)  9n€H°’°(3)  (n  -  1,2, . ); 

(b)  <pn  — >  9  as  n  — »  00 ,  uniformly  on  Ex[0,T] 

for  every  T£[0,co  ) .  Further,  suppose  ^GH1.  Then. 

9n*di|»  — >  9*di|)  as  n  — »  00,  uniformly  on  Ex[0,T]  for  every 
T€  [0,oo ) . 

Proof .  From  (a),(b)  follows  9  €C°(T3x[o,oo )) .  Choose  T€[0,oo) 
and  define  a  sequence  of  functions  fon]  on  2x(-oo,oo)  by  means  of 

=  l©n#cS|>|  ,  ©n  =  911  -  9  (n  *  1,2,  •  •  •)  •  (1.11) 

By  (a)  and  (l.ll),  ©n€  H0,®(?5) .  Moreover,  it  clearly  suffices  to 
show  that,  given  5  >  0,  there  exists  an  N(6)  such  that 

n  >  N(5 )  =>  <Cfn(x,t)  <  6  for  (x,t)  €  i5x[0,T].  (1.12) 

By  (b)  and  (l.ll),  there  exists  N(5 )  such  that  for  fixed  Tq  >  0, 

n  >  N(5)  =*  |On(x,t)|  <  5/A  for  (x,t)tRx(o,MQ],  (1.13) 

where  A  *  p^(-Tq,T)  is  the  total  variation  of  rj»  on  [-To,T],  On 
the  other  hand,  from  (l.ll)  and  a  familiar  estimate  of  Stieltjes 
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integrals  ([lO]j  p.232)  follows 


xfl(x,t)<_  max  |en|A  for  (x,  t)€ff*[0,T].  (1 .14) 

ex(o(r+-To] 

But  (1.13)  and  (l.l4)  imply  the  desired  conclusion  (1.12). 

We  have  so  far  considered  only  scalar-valued  functions 
of  position  and  time.  Functions  whose  values  are  vectors  or 
higher-order  tensors  will  consistently  be  denoted  by  underlined 
letters.  Thus,  if  the  function  v  is  defined  on  x  0,  its  value 
v(x,t)  at  position  x  and  time  t  is  a  known  tensor  for  every 
(x,t)€R  x3.  Further,  if  the  values  of  v  are  tensors  of  order 
N  >  1,  we  write  vij  ^  (N  subscripts)  for  the  components  of  v 
in  a  rectangular  cartesian  coordinate  frame  and  henceforth  adopt 
the  usual  summation  convention  for  repeated  subscripts.  We  say 
that  v€CM,N(ftx3)  or  that  if  and  only  if  the  corre¬ 

sponding  statements  hold  true  for  v^  lc»  Suppose  u  and  v 
are  tensor-valued  functions  of  the  same  order  N  >  1,  while  9  is 
scalar- valued,  and  let  u,v,9  have  the  same  domain  of  definition. 
We  then  adopt  the  notation 


u*d9  *  v 


u*dv  ~ 


Uu....k*4’’  *  vu. 


,k*avt 3 


(1.15) 


provided  the  Stieltjes  convolutions  involved  are  meaningful. 

We  may  now  turn  to  preliminaries  from  the  quasi-static 
linear  theory  of  viscoelasticity.  For  this  purpose  let  u,y, 
and  c,  with  the  components  and 


be  the  field  histories 
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of  displacement.  Infinitesimal  strain,  and  stress,  respectively . 
All  of  these  field  histories  are  to  be  regarded  as  functions  of 
position  and  time,  defined  on  3x(-cd  ,oo  ),  where  R  is  the  (regular) 
region  of  space  occupied  by  the  interior  of  the  body  in  its  unde¬ 
formed  state.  We  assume  the  body  to  be  originally  undisturbed  in 
the  sense  of  the  initial  conditions 

u^  =  =  0  on  Ex(-oo,0).  (1.16) 

Next,  we  recall  the  relevant  fundamental  system  of  field 
equations .  which  must  hold  throughout  the  space-time  domain 
Rx( -00,00) .  lhe  linearized  displacement-strain  relations  take 
the  form 

2ci 3  -  U1,J  +  u3,l  '  (1‘17) 

while  the  stress  equations  of  equilibrium  become 

"u.j +  pi  - 0  *  °ji  -  *  (1-18> 

where  stands  for  the  components  of  the  body-force  density  F. 

In  stating  the  stress-strain  relations  we  shall  make  use  of  the 
deviatoric  components  of  stress  and  strain  defined  by 

sij  =  °lj  "  3  6iJc5kk  '  eiJ  *  eiJ  “  3  6iJekk  '  ^1#19) 

in  which  denotes  the  Kronecker  delta.  Furthermore,  we  shall 
confine  our  attention  to  the  isothermal  relaxation  integral  law 
appropriate  to  linear,  homogeneous  and  isotropic,  viscoelastic 
solids.  This  constitutive  law  may  now  be  written  as 

8iJ  *  elJ*d°l  '  °kk  "  '-kk*d02 


9 


(1.20) 
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provided  G1  and  Gg,  which  are  functions  of  the  time  alone, 
designate  the  respective  relaxation  moduli  in  shear  and  isotropic 
compression. 

To  the  foregoing  initial  conditions  and  field  equations 
we  adjoin  the  boundary  conditions.  Prom  here  on  let 

st  -  °«nj  d-21' 

denote  the  components  of  the  traction  vector  S  on  a  surface  with 
the  outward  unit  normal  vector  n.  The  boundary  conditions 
governing  the  standard  mixed  problem  then  appear  as 

ui  “  ^i  on  Bi  *(“qo  >°°  ^  on  “°°  > 00 )  *  (1.22) 

in  which  and  Bg  are  complementary  subsets  of  the  boundary  B 

A  A 

of  R,  while  u  and  S  represent  prescribed  surface  displacements 
and  surface  tractions.  In  the  first  boundary-value  problem  Bg 
is  empty  and  the  displacements  are  given  on  Bx( -00 , 00 ) ;  in  the 
second  boundary-value  problem  B1  is  empty,  the  surface  tractions 
being  assigned  throughout  B*(-oo,co). 

With  a  view  toward  an  economical  statement  of  the  hypo¬ 
theses  underlying  subsequent  theorems  we  introduce 
Definition  1.3  (States,  regular  viscoelastic  states).  If  u  is  a 
vector-valued  —  and  e_,o_  are  second -order  tensor-valued  functions 
of  position  and  time  defined  on  rx(  -co , 00 ) ,  we  call  the  ordered 
array  ]  *  [u,_e,_c]  a  state  on  ftx(-oo,oo)  and  denote  by 
J(x,t)  *  Iu(x,t),e(x,t),o(x,t)]  the  value  of  J  at  (x,t). 

We  say  that  J  *  [u,_e,_o  ]  i_s  a  regular  viscoelastic  state 
on  ft x( -00,00 )  corresponding  to  the  relaxation  functions  Gv(v-1,2) 
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and  the  body-force  density  F,  and  write 


j  =  (u,e.,£  ]  €  y(KQ1,G2>Z)>  U.23) 


provided : 

(a)  u€H2'1(fd)  and  e.^P €  H0,0(e),  while  QyG  H2  with 
Qv(0)  >  0  (v  -  1,2); 

(b)  u>L*2.*£*Gv(  v  =  1*2)  on  the  Interior**  of  ft*(-oo,oo ) 
satisfy  the  field  equations  (1.17)>  (1.18),  (1.19),  (1.20); 

( c )  if  ft  is  an  exterior  region,  then,  as  x  s  |  x  I  •— >  oo , 


u(x, •)  =  OCx*1)  ,  u(x, •)  =  o(x_1) , 
l(x, •)  =  0(x"2)  ,  F(x, • )  =  0(x-3), 


1(1 .24) 


uniformly  on  [0,T]  for  every  T6  [O,oo  ) . 

If,  In  particular,  F  =  0  on  8*( -co , oo ) ,  we  write 


J  -  €7(6,0! ,02).  (1.25) 

In  (1.24),  as  In  the  sequel,  the  notion  of  "order  of 
magnitude"  Is  used  In  Its  standard  mathematical  connotation. 
Thus,  If  v  is  defined  on  ftx[0,oo),  <2  being  an  exterior  region, 
we  write  v(x*-)  =  0(*n)  as  x  — >oo,  uniformly  on  [0,T]  for  every 
T€  [0,oo)  If  and  only  If  :  given  T  >  0,  there  exist  numbers  p(T) 
and  M(T)  such  that  x  =  |x|  >  p  Implies  |v^j  ic(i,t)|  <  Mxn 

for  every  t€[0,T]. 

Conditions  (a)  and  (b),  which  are  partly  redundant  but 

mutually  consistent,  could  be  weakened  somewhat  at  the  expense 

_ — - - 

Recall  that  <2  may  be  open  or  closed. 
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of  more  elaborate  smoothness  assumptions .  Note  that  the  field 
histories  belonging  to  a  regular  viscoelastic  state  may  possess 
finite  jump-discontinuities  with  respect  to  time  at  t=0. 

Addition  of  states  and  multiplication  of  a  state  by  a 
scalar  constant  are  defined  as  follows .  Suppose  J  *  [u,e,tf]  and 
J '  =  tu'je.Sd.'  3  are  states  on  «*( -00,00 ) .  Then 


J  +  ]'  =  tu+u',e+e  ‘,d+d  »], 

\S  -  [XuAijXa]. 


(1.26) 


Further,  we  define  the  derivative  of  a  state  J  =  [u,.e,c[]  with 
respect  to  the  k-th  cartesian  coordinate  (in  a  given  coordinate 
frame)  by  means  of 


(1.27) 


provided  the  requisite  space-derivatives  ^  8111(1  djj  ^ 

exist . 


The  uniqueness  and  the  reciprocal  theorem  of  linear 
viscoelasticity  play  a  particularly  essential  part  in  the 
analysis  to  follow.  For  this  reason  we  include  here  statements 
of  both  of  these  theorems . 

Theorem  1.4  (Uniqueness  theorem).  Suppose 

&  =  tu,e.,o]ey(l?,01,02,F), 

J'  *  (u,,e',d']€‘Vn?,01,02,F). 

Further,  let 

u  =  u 1  on  B^x ( -oo  ,oo  )  ,  S  =  S'  on  B2x(-qo  ,oo  ), 


(1.28) 


(1.29) 
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where  and  are  complementary  subsets  of  B .  Then 

=  +  [w,0,0]  on  E*( -oo  ,oo  ) ,  (1.30) 

where ,  for  every  (x,t)  £  TT*(  -oo  ,  oo  ), 

w(x,t)  =  a(t)  +u(t)Ax  with  a,  uCH1,  (l.3l)5 


so  that  w  represents  an  (infinitesimal)  rigid  motion  of  the 
entire  body . 

Volterra's  [11]  proof  of  the  preceding  uniqueness  theorem 
was  spelled  out  in  detail  in  [l]  (Theorem  8.1)  with  limitation  to 
bounded  regular  regions  of  space.  The  extension  of  this  proof  to 
unbounded  (exterior)  regions  offers  no  difficulties  in  the 
presence  of  the  regularity  assumptions  (1.24)  which  imply  further 
that  §_(x,  •)  =  0(x“2)  as  x  — >  co ,  uniformly  on  every  interval  [0,Tj 
Theorem  1.5  (Reciprocal  theorem).  Let 


r  =  [u',e‘,o' ’3Or(T[,01,G2,P')  . 

Then,  for  each  t  £  ( -oo  ,  ao  ) , 

^  [ S*du 1 ] ( x, t ) dA  +  J  [F*du 1 ](x, t)dV  = 
J  [S '*du](x,t)dA  +  J  [F'*du](x,t)dV  = 

B  R  “  ”  “ 

J  [o*de  '  ](x, t)dV  =  J  [d'*de  ](x, t)dV  . 

R  R 


Ml. 32) 


(1.33) 


^  Throughout  this  paper  the  symbols  " • "  and  "a"  are  used  to 
indicate  scalar  and  vectorial  multiplication  of  vectors, 
respectively . 
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A  proof  of  this  theorem,  valid  for  bounded  regions, 
appears  in  [l]  (Theorem  7*^);  its  generalization  to  exterior 
regions  is  a  trivial  matter. 

Finally,  we  cite  a  result  that  furnishes  an  extension  to 
viscoelasticity  theory  of  the  Papkovich-Neuber  stress  functions 
in  classical  elastostatics  (see  Theorem  9-2  of  [l]). 

Theorem  1.6  (Generalized  Papkovich-Neuber  solution). 

(a)  Let  fj  be  open  and  bounded.  Let  F€H'1'il((?)  be  vector¬ 
valued  and  assume  Gv€  H2  with  Gv(0)  >  0  (v=l,2); 

(b)  Suppose  9  €  and  H^,1(R)  are  a  real -valued 

and  a  vector-valued  function,  respectively,  Buch  that 

V29  «  -  |  x-f  ,  |  f  on  fix  (-00,00  ),  (l-3^)6 

where 

f  =  F*dG^1*d(2G1+G2)_1;  (1-35) 

( c )  Define  a  state  [u,e_,dj  on  -oo ,  oo  )  b^;  means  of 

u  =  v(9+x *jJ» ) *d ( G1+2G2 )  -  ^<^(20^0^)  (1.36) 

in  con.lunctlon  with  (1.17)  and  (1 .19) ,  (l  .20) . 

Then 

4  -  [u,e,d]eY(«,01,02,F).  (1.37) 

The  completeness  of  the  foregoing  solution  to  the  field 
equations  was  established  in  Theorem  9-^  of  [l]. 


c  2 

°  7  is  the  spatial  gradient  operator  and  7  the  spatial 
Laplacian  operator. 
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2 .  Kelvin 's  problem  In  viscoelasticity  theory :  the  basic 
singular  state .  Higher-order  singular  states  . 

In  the  current  section  we  deal  first  with  the  problem  of 
a  concentrated  load  applied  at  a  point  of  a  viscoelastic  medium 
that  occupies  the  entire  space  E.  This  problem  is  the  counter¬ 
part  in  viscoelasticity  theory  of  Kelvin's  problem  in  elasto- 
statics  and  its  solution  supplies  the  fundamental  singular  solu¬ 
tion  of  the  field  equations  under  present  consideration. 

Kelvin's  solution  of  his  problem,  which  was  first 
published  without  proof  in  [12]  (1884),  was  later  deduced  by 
Kelvin  and  Tait  ([13],  p.227  et  seq.)  through  a  limit  process 
that  takes  as  its  point  of  departure  an  elastic  solid  subjected 
to  distributed  body  forces .  This  limit  definition  of  the  notion 
of  an  internal  concentrated  load  in  classical  elastostatics  is 
sketched  in  somewhat  greater  detail  by  Love  [5]  (Art.130)  and 
was  made  fully  explicit  in  [8],  which  may  serve  as  a  guide  for 
the  generalization  to  be  attempted  here.  The  theorem  to  which 
we  turn  now  is  intended  to  serve  a  dual  purpose:  first,  it  aims 
at  a  mathematically  precise  and  physically  natural  unique  charac¬ 
terization  through  a  limit  process  of  the  singular  problem  at 
hand;  second,  it  furnishes  the  explicit  solution  to  the  problem 
thus  formulated.  To  facilitate  the  statement  of  this  theorem  we 
first  agree  that  Qp(x°)  denotes  an  open  sphere  with  the  radius  p 
centered  at  x°  and  adopt 

Definition  2.1  (Limit  definition  of  a  concentrated  load).  We  say 
that  [if1]  is_  a  sequence  of  body-force  distributions  that  tends  to 
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a  concentrated  load  L  applied  at  x°  if  x°€  E  and  L  e-H^  Is 
vector-valued,  while  {Ff1}  has  the  following  properties : 

(a)  for  every  n  (n=l,2,  •  •  • )  if1  is  a  vector-valued  function 

with 

FfeH2'1^)  ,  Pri  =  0  on  (E-Qn)v(-oo  ,oo ),  (2.1) 

where  {Qn}  .is  a  sequence  of  "load  regions"  characterized  by 

Qn  =  Q  n(x°)  ,  6n  — >  0  as  n  — >  oo  ;  (2.2) 

5  n 

(b)  J  f(x,  -)dV  ->  L  as  n  — >  oo  ,  uniformly  on  [0,T] 

Qn 

for  every  T€[0,oo); 

(c)  the  sequence  of  functions  defined  by 

4>n  -  |  l|?(x,-)|dV  on  (-00,00)  (n=l,2,-..),  (2.3) 

Qn 

Is  uniformly  bounded  on  [0,T]  for  every  T€[0,oo). 

We  are  now  in  a  position  to  turn  to 
Theorem  2.1  (Viscoelastic  Kelvin-state).  Let  (Pn|  be  a  sequence 
of  body-force  distributions  that  tendB  to  a  concentrated  load  L 
applied  at  3 j° .  Let  Qv€  H2  with  Gy(0)  >  0  ( v=l,2) .  Then: 

(a)  there  exists  a  unique  sequence  of  states  £|n}  SHSU 

that 

4n  =  kn,.1n,dn]€^(E,G1,Q2,Fn)  (n=l,2,  •  •  • ) ;  (2.4) 

(b)  4  n  converges  to  a  limit  state  |  as  n  — >  00 ,  uniformly 
on  ?5x(-oo,T]  for  every  bounded  ft  such  that  x°^  and  every 

T  €  ( -00  ,00  ) ; 
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(c)  the  limit-state  J  *  [u,_e, d ]  is.  Independent  of  the 
particular  choice  of  the  sequence  (FnJ  and  —  in  the  sense  of 
Theorem  1.6  —  Is,  generated  by  the  stress  functions  nj»  defined 
through 


<p(x,t)  =  0  ,  ±(x,t)  = 


f(t) 


8x|x-x°| 


(2.5) 


for  every  x  /  x°  and  every  t  £ ( -oc , oo ) ,  where 

f  =  L*d0"1#d(201+G2)"1  on  (-oo,co).  (2.6) 

We  call  ^  the  (viscoelastic)  Kelvin-state  corresponding 
to  a  concentrated  load  L  applied  at  x°  and  to  the  relaxation 
functions  G-^Gg. 

Proof.  Note  that  the  existence  of  the  Stieltjes  Inverses  G"1, 
(2G1+G2)"1  is  assured  by  Theorem  1.2  and  the  present  hypotheses 
on  0v(v=l#2).  Define  two  sequences  of  functions  {tpn}  and 
by  setting,  for  every  n(n*l,2,-**)  and  all  (x,t) €  Ex( -go  ,oo ), 

L'£n(i>t) 

Qn 

n/  x  .if  In(L,t) 

1  (x*t)  =  w  J  —  dV, 

Qn 


n,  If  L'£n(i>t)  _ 

<P  (i,t)  =  Bx  i  n  |  x'-Ll —  A\  * 

on  1  —  3- 


k  • 


(2.7) 


where 


fn  -  Fn*dG^1*d(2G1+G2)-1  on  Ex(-oo,oo  ) .  (2.8) 


A  subscript  attached  to  an  "element  of  volume"  or  an  "element 
of  area"  in  a  volume  or  surface  integral  indicates  the 
appropriate  space  variable  of  integration. 
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Thus  <pn(‘,t)  and  tjin(*,t)  are  Newtonian  potentials  of  mass  distri¬ 
butions  over  If1.  It  follows  from  (a)  in  Definition  2.1, 

Theorem  1.1,  and  a  trivial  extension  of  Lemma  9*1  in  [1]  that 
<pn,  ]J)neH3,1(E)  and  that 

V2<pn  =  -  -  x-fn  ,  V24[n  -  |  fn  on  E*( -00,00 ).  (2.9) 

Consequently  and  because  of  the  behavior  of  the  integrals  in 
(2.7)  as  x  — ■»  00,  the  functions  <pn  and  ]j>n,  v/hen  used  as  stress 
functions  in  conjunction  with  Theorem  1.6,  generate  a  state 
=  [jinjen,fln]  that  meets  (2.4).  The  uniqueness  of  $n  is 
immediate  from  Theorem  1.4.  Thus  part  (a)  of  the  theorem  under 
consideration  has  been  confirmed. 

Let  J  =  [u,je,^j]  be  the  state  generated  by  the  stress 
functions  in  (2.5).  Choose  ft  such  that  x°^  ft,  choose 
T£(0,oo),  and  hold  K'<(0,T]  fixed  for  the  remainder  of  the 
argument.  In  order  to  establish  parts  (b)  and  (c)  of  the 

g 

theorem  It  suffices  to  show  that,  as  n  — >00, 

un  — >u  ,  _en  — > _e  ,  011  — >  c_  uniformly  on  ftx[0,T].  (2.10) 

By  virtue  of  Theorems  1.6,  1.3*  however,  (2.10)  follows  if  we 
show  that  for  a  fixed  choice  of  the  coordinate  frame,  as  n  — >  00 
and  uniformly  on  Rx[0,Tj, 

0  Observe  that  both  i  and  vanish  identically  on  Wx(-oo,0). 
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9n  -»  9  ,  ->  V  i 

^in  ->\ji  ,  ^  A 


_n  v  _ 

9,ij  ->9.u  ■ 


<-  (2.11) 


Since  all  of  (2.11)  may  be  established  by  strictly 
analogous  means,  we  shall  demonstrate  merely  that 

—>  ^  as  n  — >  oo  ,  uniformly  on  ft*[0,T).  (2.12) 

To  this  end  we  infer  from  (2.5)>(2.7)  that  for  (x,t)  €  ft*[0,T], 

8nfin(x,t)  -  £(x,t)]  =  lj(x,t)  +  l£(x,t),  (2.13) 


provided 


-  I  i-~7  -  -i—)fn(i,t)dV 
Qn  |x-x|  |x-£l  i 

l£(x,t)  =  -  -1  -  ■  (f(t)  -J  fn(£,t)dV]  . 

|x-xl  Qn 


y  (2.14) 


Accordingly,  it  is  sufficient  to  show  that 

— 1J— 2  0  as  n  00  •  unifonnly  on  «x[0,T].  (2.15) 

In  view  of  (2.8),  the  first  of  (2.14)  becomes 

-  I  (— i—  -  -i_][F",<i0](S.,t)av,  (2.16) 

QTl  I  X-X  I  |x-£| 

for  all  (x, t) £ftx[0,T],  if  one  sets 

0  =  0‘1*d (20^02 r1  on  (-00,00).  (2.17) 

Now,  by  the  hypothesis  underlying  the  present  theorem  and  (a)  in 
Definition  2.1,  there  is  an  N  >  0  such  that  ‘SnQ*1  is  empty  when 
n  >  N.  Consequently,  for  n  >  N  and  (x,t)  €  fgx[0,T),  the  volume 
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integral  (2.16)  is  proper.  This  entitles  us  to  reverse  the  order 
of  the  process  of  space-integration  and  convolution  in  (2 .16) 
when  n  >  N,  as  is  readily  seen  with  the  aid  of  Theorem  1.2,  (a) 
and  (g)  in  Theorem  1.1,  and  the  available  regularity  of  the 
functions  Pn  and  Gv(v=i,2).  Hence 

ij  =  In*dG  on  13x[o,T]  (n  >  N),  (2.18) 

where 

In(x, t)  =  f  [— i-s - —  jF^tJdV,  (2.19) 

li-x  |  |x~Ll  L 

and  clearly  In  is  continuous  on  ftx[0,Tj.  From  (2.l8)  follows  the 
estimate,  valid  for  n  >  N  and  fixed  Tq  >  0, 

|I?|  <_  max  |In|  nrr(-T  ,T)  on  «x[0,T],  (2.20) 

1  ftx[0,T+To]  u  ° 

where  Wq(-Tq,T)  is  the  total  variation  of  G  on  [-Tq,T].  On  the 
other  hand,  one  draws  from  (2.19)  that 


I  In(x,t)| 


< 


max 


"  ~  r|  I^I.tJIdV 

l*-il  I  Qn 


(2.21) 


for  every  (x»t)  £  f£x[o,T4-To] .  But,  in  view  of  (c)  in  Defini¬ 
tion  2.1,  the  integral  in  (2.21)  is  bounded  uniformly  for  all  n 
and  all  t€[0,T+T  ],  whereas  the  coefficient  of  this  integral 
tends  to  zero  uniformly  as  n  — >  co  for  all  x£fc  because  of  (a) 
in  Definition  2.1.  Consequently  I11  — >  0,  as  n  — »co,  uniformly 
on  Rx[o,Th-Tq]  and  hence  (2.20)  implies  the  statement  concerning 
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Integral  (2.16)  Is  proper.  This  entitles  us  to  reverse  the  order 
of  the  process  of  space-integration  and  convolution  in  (2.16) 
when  n  >  N,  as  is  readily  seen  with  the  aid  of  Theorem  1.2,  (a) 
and  (g)  in  Theorem  1.1,  and  the  available  regularity  of  the 
functions  Fn  and  0v(v=l,2).  Hence 

l£  =  ,In*dG  on  T3x[o,T]  (n  >  N),  (2.18) 

where 

In(x,t)  =  f  [_L_ - JPn(L,t)dVE  (2.19) 

Qn  |x-x°|  |x-Ll  1 

and  clearly  In  is  continuous  on  ftx[0,T].  Prom  (2.18)  follows  the 

estimate,  valid  for  n  >  N  and  fixed  T  >  0, 

o 

|I?|  <  max  |In|  7Tr(-T  T)  on  ®x[o,T],  (2. 20) 

1  «x[0,T+Tq]  u  0 

where  Vq(-Tq,T)  is  the  total  variation  of  G  on  [-Tq,T].  On  the 
other  hand,  one  draws  from  (2.19)  that 

|Xn(x,t)l  <  max  ,  1  -  — M  f  |Fn(i,t)|dV  (2.21) 

“ie5nls-2°l  l2-S.llJcn 

for  every  (x,t)  €Rx[0,T+Tq]  .  But,  in  view  of  (c)  in  Defini¬ 
tion  2.1,  the  integral  in  (2.21)  is  bounded  uniformly  for  all  n 
and  all  t€[0,TM-T  ],  whereas  the  coefficient  of  this  integral 
tends  to  zero  uniformly  as  n  — >  co  for  all  x£  fc  because  of  (a) 
in  Definition  2.1.  Consequently  I11  — >  0,  as  n  — >  co ,  uniformly 
on  Rx[o,TH-To)  and  hence  (2.20)  implies  the  statement  concerning 
ij  in  (2.15). 
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By  virtue  of  (2. 6), (2. 8),  and  (2.17)*  the  second  of  (2.14), 
after  a  permissible  reversal  of  the  processes  of  space-integration 
and  convolution,  yields 

Io(x,t)  =  |-^-0--r{L  -  j  jf(I,-)dV}*dO](3C,t)  (2.22) 

|x-2i  I 

for  all  (x,t) €  ftx[0,T] .  The  assertion  concerning  l£  in  (2.15) 
now  follows  from  (b)  in  Definition  2.1  together  with  Theorem  1.3. 
This  completes  the  proof . 

Conditions  (a)  and  (c)  in  Definition  2.1  trivially  imply 

(c')J  (x-x°)AFn(x,  -)dV  -*0  as  n  — >  oo,  uniformly  on  [0,T] 

Gn 

for  every  T€[0,oo).  It  is  natural  to  ask  whether  the  seemingly 
artificial  requirement  (c)  may  be  replaced  by  (c!)  without 
impairing  the  conclusion  in  Theorem  2.1.  That  this  is  not 
possible  is  clear  from  Theorem  4.3  in  [8],  on  passing  to  the 
special  case  of  the  elastic  solid.  Indeed,  If  (c)  is  replaced 
by  (c1)  In  Definition  2.1,  one  can  construct  a  sequence  {if1} 
of  body-force  distributions  that  tends  to  a  load  L  **  0  at  x° 
and  yet  generates  —  in  the  sense  of  the  limit  process  underlying 
Theorem  1.2  —  a  limit  state  other  than  the  null-state.  Require¬ 
ment  (c)  does  become  superfluous,  however,  if  (if1}  Is  restricted 
to  body  forces  that  are,  at  every  instant,  parallel  and  of  the 
same  sense.  In  this  special  case  one  has 

Fn(x, t)  =  I Fn(x,t)j k(t)  for  all  (x,t) €  Qnx[0,oo  ),  (2.23) 

k(t)  being  a  unit-vector,  and  (a),(b)  are  readily  found  to  imply 
(c),  in  Definition  2.1. 
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Prom  here  on  let  denote  a  unit -vector  in  the  xa- 
direction.  If  L  =  hea  In  Theorem  2.1,  we  call  the  limit  state 
J  the  normalized  (viscoelastic)  Kelvin-state  (corresponding  to 
a  load  in  the  xa-direction  applied  at  x°  and  to  the  relaxation 
functions  G^G^).  Moreover,  we  denote  the  value  of  this  state 
at  (x,t)  by 


with  the  understanding  that  all  superscripts  not  otherwise 
specified  henceforth  have  the  range  of  the  integers  (1,2,3)- 
Evidently,  for  every  constant  vector  a, 

<T(x,t;x°+a)  »  Ja(x-a,t;x°),  (2.25) 

whence  in  particular 

vfa(x,t;x°)  =  /°(x-x°,t;0) .  (2.26) 


We  now  record  the  cartesian  components  of  displacement 
and  stress  belonging  to  fa(x,t;0) .  These  components  are  easily 
obtained  on  setting  L  =  h^  and  x°  *  0  in  (2.5),  by  recourse  to 
(1.36),  (1.17),  (1.19),  and  (1.20).  In  this  manner  one  arrives  at 

X_x,  XaX, 

X-  -  -  x‘ 

«  -A(  — X-fi[2h(t)  -  3Qp(t)] 

SltX'3  d 


C»al  +  52^1  +  3«j(t)[6ol  -  ] 


M2. 27) 


where  x  =  |x|  =  »  as  before,  and  ar®  auxiliary 

response  functions  defined  by 
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J1  =  °i1  '  Q1  =  (SQ^)"1  >  Qg  =  VdGl  on  (-oo,oo).(2.28)9 
In  the  particular  case  of  an  elastic  solid  one  has 


provided  p.  and  x  stand  for  the  shear  modulus  and  the  bulk 
modulus  of  the  material,  respectively.  In  this  instance  (2.27), 
for  all  x  7  0  and  0  t  <00,  reduce  to  the  analogous  formulas 
appropriate  to  the  normalized  elastostatlc  Kelvin-state  (see, 
for  example,  equations  (5.4)  in  [8]). 

We  take  up  next  a  theorem  that  summarizes  certain 
relevant  properties  of  J°(x,t;0).  For  the  sake  of  convenience 
in  stating  this  and  subsequent  results  we  introduce  the  symbol 

I 

E  for  the  open  region  consisting  of  all  points  of  the  Euclidean 
x 

space  E  with  the  exception  of  the  point  x°.  Further,  we  shall 

I  O 

simply  write  E’  in  place  of  E  when  x  =  0,  so  that  E*  denotes 

x 

the  complement  of  the  origin  with  respect  to  the  entire  space  E. 
Theorem  2.2  (Properties  of  the  Kelvin- state ) .  The  normalized 
Kelvin-state  Ja(x,t;0)  has  the  properties : 

(a)  ua,ia,ia£H®,1(E')  and  Ja  -  [u°,ia,Ao]€*Y(E,,01,(^]; 

» 

(b)  .  Sa(x,  •  ;0)dA  -  e^  on  (-00,00), 

where  n  is  any  closed  regular  surface  surrounding  the  origin  and 

9 


Note  that  J,  is  the  creep  compliance  in  shear  (see  Theorem  3.3 
in  [1]) . 
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Sa  is  the  traction  vector  of  on  that  side  of  n  which  faces  the 
origin: 

(c)  |  xa  Sa(x#-jO)dA  =  0  on  ( -oo  ,  oo  ) ; 

n 

(d)  ua(x,  *;0)  =  0(x-1)  ,  o.a(xj  *  ;0)  =  0(x“2)  as  x  — »  0, 
uniformly  on  [0,T]  for  every  T  € [0,oo  ) . 

Proof.  Property  (a)  is  Immediate  from  the  specific  form  of  the 
stress  functions  generating  exhibited  in  (2.5),  and  from 

Theorem  1.6.  To  confirm  (b),  use  the  second  of  (2.27)  and  note 
that  because  of  (a)  the  surface  II  may  be  restricted  to  be  a 
sphere  centered  at  the  origin.  Property  (d),  which  characterizes 
the  order  of  the  displacement  and  stress  singularities  at  the 
load-point  x°  *  0,  follows  at  once  from  (2.27),  whereas  (c)  is 
readily  found  to  be  implied  by  (a)  and  (d). 

In  the  treatise  literature  on  elasticity  theory  the 
original  Kelvin-problem  is  frequently  formulated  on  the  basis 
of  the  elastostatic  counterpart  of  (a),(b),(c)  in  Theorem  2.2. 

As  was  pointed  out  in  [8],  such  a  direct  formulation  of  the 
problem  fails  to  determine  its  solution  uniquely.  Similarly, 
(a),(b),(c)  do  not  furnish  a  unique  characterization  of  the 
viscoelastic  Kelvin-state  J°(x,t;0)  since  —  as  we  shall  have 
occasion  to  see  shortly  —  there  exist  viscoelastic  states 
(distinct  from  the  null-state)  that  are  regular  on  E'x(-co,a>) 
and  possess  self -equilibrated  singularities  at  the  origin.  In 
this  connection  we  mention  further  that  the  viscoelastic  Kelvin- 
state  i°(x,tjO),  whose  displacements  and  stresses  are  exhibited 
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in  (2.27),  may  alternatively  be  obtained  directly  from  the  corre¬ 
sponding  elastostatic  state  through  a  purely  formal  application 
of  the  well-known  correspondence  principle  that  links  the  linear 
theories  of  elasticity  and  viscoelasticity.10  However,  this 
elementary  method  for  deducing  the  viscoelastic  Kelvin-state 
does  not  assure  the  truth  of  Theorem  2.1,  from  which  the  singular 
state  in  question  derives  its  intrinsic  physical  significance. 

We  turn  now  to  a  discussion  of  the  higher-order  singular 
viscoelastic  states  that  may  be  generated  through  a  single  space- 
differentiation  of  the  Kelvin-state.  Thus  let  Ja(x,t;x°)  once 
again  be  the  normalized  viscoelastic  Kelvin-state  introduced 
previously  and  define  a  set  of  nine  viscoelastic  doublet -states 

f (x,t ;x°)  =  [u0^(x,tiX0),£a^(x,t;x0),d°P(x,t,*x0)  ]  (2.30) 

by  means  of 

JaP(x,tjx°)  =  Ja(xit;x°)  for  (x,t)€  e'ox(-oo,co  )  .(2.31)11 
P  X 

A  physical  interpretation  of  is  easily  established.  Indeed, 
(2.31)  implies 

|a^(x>t;x°)  =  lim  [-gt  Ja(x+©eR,t  jx°)  -  J a(x,t ;x°)  3] ,  (2 .32) 

10  For  the  analogous  treatment  of  the  problem  of  the  half -space 
under  a  concentrated  surface  load  perpendicular  to  the 
boundary,  see  Lee  [14]. 

11  Recall  (1.27). 
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which,  because  of  (2.25),  is  equivalent  to 

=  llm  lsUa(^t;x°“©eQ)  ~  |a(x>t;x°)  ]}  .(2.33) 
0  ©  _»  0  P 

Prom  (2.3l)  and  (2.26)  follows 

fP(x,t;x0)  =  JGP(x-x°,t;0).  (2.34) 

We  list  next  the  cartesian  components  of  displacement  and  stress 
belonging  to  Ja^(x,tjO),  which  may  be  computed  from  (2.27)  with 
the  aid  of  the  defining  relation  (2.31). 

u“P(x,t;0)  =  {[2J1(t)+3Q1(t))6laXp 

+  l2Jl<t>-3Vt>H^££  -  ^  -  ^xj}  , 

X 

OlPj(x,tiO)  -  ^-5  tt2h(t)-3Qg(t)][-X^X.l  -  6plxax3  -  (2.35) 

-  ‘pjVl'Wj1  +  Q2(t,[‘>ai(3xpxJ-,>Mx2> 

+  ^ojOXiXp-^x2)  -  '■1J(3xaXp-6(1px:3)]j  . 

Here  J^,Q^,Qg  are  again  the  response  functions  given  by  (2.28). 

The  subsequent  theorem  is  a  trivial  consequence  of  (2.31), 
Theorem  2.2,  and  the  explicit  formulas  (2.35). 

Theorem  2.3  (Properties  of  the  doublet -states ) .  The  doublet -state 
J“P(x,t;0)  ha3  the  properties : 

(a)  aaP>£aP  a«P€H°o.l(E.)  and 

!a?  *  [u°P,t°P,0aP]6  'Y(E',Q1,a2); 

* 

(b)  Sa^(x,  •  ;0)dA  =  0  on  (-00,00), 
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where  n  is.  any  closed  regular  surface  surrounding  the  origin  and 
S°P  is  the  traction  vector  of  JaP  on  that  side  of  n  which  faces 
the  origin? 

(c)  J  xa  SaP(x,  •  ;0)dA  =  6  Q  eh  on  (-00,00), 

jj  -  _  “  apiTY  — 

where  designates  the  components  of  the  usual  alternator; 

(d)  u°P(x, -;0)  =  0(x~2),£°P(x,*j0)  =  0(x'3)  as  x  — >  0, 

unlfomly  on  [0,T]  for  every  T€  [0,oo  ) . 

Properties  (a),(b),(c)  were  to  be  anticipated  Intuitively 
because  of  the  physical  meaning  attached  to  the  doublet-states 
by  (2.33)j  and  in  view  of  (a), (b), (c)  of  Theorem  2.2.  As  Is 
apparent  from  (b)  and  (c)  in  the  present  theorem,  the  stress 
singularity  of  Ja^(x,UjO)  at  the  origin  13  statically  equivalent 
to  a  couple  or  to  null  depending  on  whether  a  /  p  or  a  =  p.  For 
this  reason  we  call  ^a^(x,t;x°)  a  doublet -state  with  moment  or  a 
doublet -state  without  moment,  according  as  a  /  p  or  a  =  p. 

As  In  elasticity  theory,  it  is  expedient  to  introduce 
two  particular  linear  combinations  of  doublet -states .  Thus  we 
designate  by 

J°(x,tjx°)  =  [u° ( x , t ;x° ) ,e. 0 ( x> t ;  x° ) ,g_° ( x, t ; x° )  ]  (2.36) 

the  state  defined  through 

l°(x,tjx°)  =  |aa(x,t;x°)  (2.37)12 

The  summation  convention  is  henceforth  understood  to  apply  also 
to  repeated  superscripts  and  to  repeated  indices  that  appear 
once  as  a  subscript  and  once  as  a  superscript,  provided  these 
indices  have  the  range  (1,2,3). 


12 
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and  call  It  the  state  appropriate  to  a  (viscoelastic)  center  of 
compression  at  x°.  Further,  we  denote  by 


=  lSa(x,t;2S°)iIa(x>tjx0)15a(x,t;x0)  ]  (2.38) 


the  state  defined  through 


apY' 


(2.39) 


which  we  address  as  the  state  appropriate  to  a  (viscoelastic) 
center  of  rotation  in  the  xa-direction  at  x° . 


From  (2.37),  (2.39),  (2.35)  follow 
u?  (x,tjO)  =  x.Q,  (t)  , 


“l  *  ^3  6clfVl(t)  ’ 


(2.40) 


>(2 .41 ) 


Ihe  implications  of  Theorem  2.3,  as  far  as  the  properties  of  the 
states  appropriate  to  a  center  of  compression  and  a  center  of 
rotation  are  concerned,  are  immediate.  In  particular  J°(x,t;0) 
has  a  3elf-equilibrated  singularity  at  the  origin,  whereas  the 
singularity  inherent  in  J^^tjO)  is  statically  equivalent  to  a 
couple  whose  axis  is  the  xa-axis.  Specifically, 
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Through  successive  space-differentiations  of  the  Kelvin- 
state  /a(x,t;x°)  one  may  evidently  generate  an  infinite  aggregate 

i  . 

of  viscoelastic  states  that  are  regular  on  E  x(-co,oo)  and 

x° 

possess  singularities  at  x°  of  progressively  higher  order. 
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3.  Qreen's  states.  Integral  representations  for  the  solution 
of  the  fundamental  boundary-value  problems . 

As  a  prerequisite  for  the  treatment  of  the  main  subject  of 
this  section  we  require  two  lemmas  concerning  certain  integral 
properties  of  the  viscoelastic  Kelvin-  and  doublet -states  dlB- 
oussed  in  the  preceding  section.  To  shorten  the  formulation  of 
these  lemmas,  and  for  future  convenience,  we  make  the  agreement 
that  Sp(x°)  henceforth  denotes  a  spherical  surface  of  radius  p, 
centered  at  the  point  with  the  position  vector  x°.  Further, 
Qp(x°),  as  before,  denotes  a  epherical  neighborhood  of  x°,  with 
the  radius  p.  Finally,  we  shall  write  Zp  and  Qp  in  place  of 
£p(x°)  and  Qp(x°),  when  x°  »  0. 

Lemma  3.1.  Let  e  be  a  neighborhood  of  a  point  x° €  E .  Let 

[u,e,o]G'Y(«,01,02,F) 

and  suppose  J®(x,t;x0)  is  a  normalized  viscoelastic  Kelvin-state . 
Then  for  fixed  x°  and  each  t  6-  ( -co  ,  co  ) : 


(a) 

11m 

J  [S*dua](x,t)dA  *  0, 

(3.1) 

p  — >  0 

2p(x°) 

(b) 

11m 

f  [Sa*du](x,t)dA  *  ua(x°,t). 

(3.2) 

p  — >  0 

2p(x°) 

provided  S  and  Sa  are  the  respective  traction  vectors  of  /  and 
|a  on  that  side  of  Sp(x°)  which  faces  x°. 

Proof .  Without  loss  in  generality  assume  x°  ■  0.  Also,  since 
|(x,t)  and  fa(x,t;0)  vaniBh  for  all  (x,t)6 ftx(-oo  ,0),  it 
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suffices  to  consider  t  ^  0.  Thus  choose  t€  (0,oo)  and  hold  it 

fixed  throughout  the  remainder  of  the  argument.  Further,  let 

p  >  0  be  such  that  ZL  C  fc  . 

°  Ho 

With  a  view  toward  proving  (a),  define  I°(p)  for  all 

/ 


P£(0,pQJ  by  means  of 


lT(p)  *  r.  [S*dua ](x,t)dA, 

1  EP 

whence  and  fi'om  (b)  in  Theorem  1.1, 

l*l(p)l  <  |  |  [\/I*dS](x,t)|dA. 
EP 

Now  note  that 


(3.3) 


(3.4) 


’KH1  =^irr(-oo,t)  <  U(o)|  +  t  max  U|  ,  (3.5) 

w  [0,tj 

where  -oo, t)  is  the  total  variation  of  t|>  on  (-co,t].  Prom 

(3-4),  (3.5)*  in  view  of  (1.15),  (1.8),  and  the  usual  estimates 
of  Stleltjes  and  Rlemann  Integrals,  follows 


1 1? ( P ) I  <  I  max  |u®| ][max|s( *,0)|  +  t  max  |s|]4xp2.  (3.6) 

1  2px[o,tJ  2p  “  2px[0,t)~ 

The  leading  term  within  brackets  is  0(p_1)  as  p  — »  0  because  of 
(d)  in  Theorem  2.2,  whereas  the  second  term  within  brackets  is 
uniformly  bounded  fpr  all  p€(0, pQ]  by  virtue  of  (1.21)  and  the 
regularity  of  x. 1  on  3  *[0,tj.  Hence  (3.6)  implies  (3.1). 

J  H°  n 

Consider  next  part  (b)  and  define  Ig(p)  for  all  p  €(0,pQ] 
through 

I !(p)  -  /  [Sa#du)(x,t)dA. 

EP 


(3.7) 
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Setting 


v(x,-t)  -  u(x,t)  -  u(0,t) 


(3.8) 


for  all  (x,'t)Gfik(-oo ,oo )  and  invoicing  once  more  (b)  in 
ftieorem  1.1,  one  has 

llo(p)  -  ua(0,t)|  <  I  J  [v#dSa](x,t)dA| 

2P 

+  |[u(0,-)*dj  Sa(x,sO)dA](t)  -  ua(0,t)|  .  (3.9) 

2P 

The  second  term  in  the  right-hand  member  of  (3.9)  vanishes  by 
virtue  of  (b)  of  Theorem  2.2  and  (f)  of  Theorem  1.1.  To 
estimate  the  first  term  proceed  as  in  part  (a) .  In  this  manner 
one  is  led  to 

ll£(p)  -  ua(0,t)l  < 

[  max  | v|  Umax  |S°(*,0;0)|  +  t  max  |Sa|]4*p2.  (3.10) 

Zpx[0,t]  Zp  Zp*[0,t]  ” 

The  leading  term  within  brackets  is  o(l)  as  p  — >  0  because  of 

(3.8)  and  the  continuity  of  u  on  B  x[0,t).  On  the  other  hand, 

PD 

— p 

the  second  term  within  brackets  iB  0(p  )  in  view  of  (l.2l)  and 

(2.27).  Thus  (3*2)  follows  and  the  proof  of  Lemma  3.1  is 
complete. 

Lemma  3.2.  Let  ft  be  a  neighborhood  of  a  point  x°e  E.  Let 

[u,e,o)fe'V(^01,02,P) 

and  suppose  Ja^(x,t;x°)  is,  a  viscoelastic  doublet -state  corre¬ 
sponding  to  the  relaxation  functions  C^, Qg .  Then  for  fixed  x° 
and  each  t  G  (  -od  ,  00  ) : 
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(a)  lim  J  [S#duaP)(x,t)dA  -  |  eap(x°,t)  +  Xao(x°,t), 


P — *  0  Sp  (x°) 


(3.11) 


(b)  lim  J  [SaP*du](x,t )dA  = 
p — >  0  Sp  (x°) 

^p>a(x°,t)  -  |  eap(x°>t)  +  Xap(x°,t),  (3.12) 


where 


\,p<2°.t)  -  5  [iapeYY(2°.t)  +  2t<cap-2!'0?eYY)*IiQ2)(-0’t)}  * 

( 3 • 13 ) 


while  S  and  SaP  are  the  respective  traction  vectors  of  $  and 
|aP  on  that  3lde  of  2p (x°)  which  faces  x° • 


Proof .  As  in  the  proof  of  Lemma  3-1j  assume  x°  *  0,  hold 

t£  [0,oo)  fixed,  and  let  p  >  0  be  such  that  Z  C  fi  , 

°  aB  P° 

Consider  first  part  (a)  and  define  I1P(p)  for  all 
p€(0,  p  ]  through 


I?^(p)  *  J*  [S*du0^ ](x,t)dA . 
J-  ? 


(3.14) 


In  view  of  the  present  hypotheses,  one  draws  from  (3.14)  in 
conjunction  with  (2.35)>  (1.21),  the  divergence  theorem. 
Theorem  1.1,  and  (1.18),  after  a  straightforward  computation. 


that 

if  (p)  -  i“p(p) +  i“p(p),  (3.15) 

where 

!®P(p)  =  _1  J  ([d  B*d(2J,+3Q1)](x,t)  + 

d  8np 0  Qp  H 
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rf(p)  =  -  I  {[P^dtaJi+SQJKx^tJxg  + 

3  8*p3  np  1  1  P 

+  [F1«d(2J1-3Q1)  ](x,t)  [— -  6apX1-6pixa]}dV.  (3.17) 

By  Definition  1.3  and  (a)  in  Theorem  1.1,  the  convolutions 
entering  the  integrand  in  (3.16)  are  (for  fixed  t)  continuous 
on  Qp,  whereas  the  integrand  in  (3.17)  is  o(l)  as  x  — »  0  so  that 
I°P(p)  =  o(l)  as  p  — »  0.  Consequently  (3.15),  (3.16),  (3.17) 
furnish 

1°|!(p)"i^a  (toafi*<i<2Ji+3V)<0-t> 

+  Iolj*d<2J1-3ei)l(0,t)[^(613xaxp«I1Jx1Xp-H.pjxixa) 

-  +  o(l)  as  p-»0.  (3.18) 

Now  carry  out  the  space  integration  in  (3.18)  and  use  (d)  of 
Theorem  1.1,  as  well  as  the  second  of  (l.l8),  to  obtain 

lf(p)  -  -J5  (6[Pap.<i(J1+(i1))(0,t) 

•  ^pf<5ii*d(2Ji-3Ql)J(°'t)}+  o(l)  as  p  >  0.  (3.19) 

Conclusion  (a)  follows  from  (3.19),  (3.1*0  after  an  elementary 
computation  involving  the  use  of  Theorem  1.1,  Theorem  1.2, 
equations  (1 .19) , (1 .20), (2 .28) ,  and  the  definition  of  X^given 
by  (3.13). 

Turn  to  part  (b)  and  define  I^(p)  for  all  p€(0,PQ] 


through 
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I?P(P)  =  l  ISap#du](x,t)dA.  (3.20) 

g  p 

The  required  limit  of  I^p),  as  p  — >  0,  may  be  confirmed  by 
an  argument  which  is  quite  similar  to  that  employed  in  establish¬ 
ing  part  (a) .  In  this  connection  one  needs  to  make  use  of  the 
fact  —  implied  by  (a)  in  Definition  1.3  —  that,  for  each 
(x,t)€  Qpx[0,oo  ), 

u(x,t)  =  u(0,t)  +  u  .(0,t)x.  +  U(x,t)  (3*21) 

)  J  J 

with  U€C2jl(Qpy[0,oo))  and 

U(x,  * )  =  0(x2)  as  x -»  0,  (3.22) 

uniformly  on  [0,t].  Further  details  of  the  proof  may  safely  be 
omitted . 

Lemma  3.2,  by  virtue  of  (2.39)*  at  once  yields 
Corollary  3.1.  Let  x°,G,  and  j  meet  the  same  hypotheses  as  in 
Lemma  3.2.  Suppose  J^^tjx0)  i_s  the  state  appropriate  to  a 
center  of  rotation .  Then  for  fixed  x°  and  each  t  € ( -oo , oo )  : 

(a)  lim  J  (S^dx^KxjtJdA  =  0, 

P->0  Zp(x°)  " 

(b)  lim  J  ['Sa#du)(x,t)dA  *  wrt(x°,t), 

P->0  2p(xO)  -  ~  “ 

where  w  is  the  rotation  field  history  of  J  ,  1 .e. 

u  =  -|  VA  u  on  Kx(-oo  ,oo  ),  (3.23) 

while  S  and  1°  are  the  respective  traction  vectors  of  /  and 
Ja  .22  that  side  of  Zp(x°)  which  faces  x°. 
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It  is  worth  mentioning  that  conclusions  (a),(b)  in 
Lemma  3*1  and  Corollary  3.1  do  not  involve  the  relaxation  func¬ 
tions  of  either  the  singular  or  the  regular  state  under  consider¬ 
ation.  Also,  the  conclusions  in  Lemma  3.1  and  Corollary  3.1 
hold  true  even  if  the  relaxation  functions  of  the  singular  state 
are  distinct  from  those  belonging  to  the  regular  state. 

We  are  now  in  a  position  to  deduce  the  integral  repre¬ 
sentations  which  constitute  the  main  objective  of  this  section. 
With  a  view  toward  a  representation  theorem  applicable  to  the 
first  fundamental  boundary-value  problem,  in  which  the  surface 
displacements  are  prescribed  over  the  entire  boundary  for  all 
time,  we  introduce 

Definition  3.1  (Green’s  states  of  the  first  kind).  We  call 

fa(£»t;x)  =  [ua(^,t;x),la(^,t;x),3a(Jl,tjx)  ], 

Jap(]L*t;x)  =  [u°P(^,t  jx),i°^(i,t;x),3aP(^,t;x) ) 


the  Green’s  states  of  the  first  kind  for  a  (regular)  region  R 
and  relaxation  functions  G^Og  IX  and  only  if:  for  all 
(j^,t;x)  €"RX( -oo,oo )xR  with  £  ft  x, 


/a^(i»t;x)  = 


J°(£*t;x)  =  J°(i,t;x)  +J°(£,t;x), 

-  !UaP(£,t;x)  +  )P<X(I,t;x)]  +/aP(I,t;x), 


M3.24) 


where  /a(£,t;x)  and  7a^(£,t;x)  are  the  normalized  Kelvin-state 
and  the  doublet-state,  corresponding  to  01 , Qg ;  further 


j6z( 25)/21 


38 


|a(L>  tjx)  =  [ua (£,  t  ;x)  ,6° (£,  t  jx)  ,3° (£,  t  ;x)  ] , 

for  each  x£R,  are  states  such  that 

(a)  i“P(-.-J5)6mo1.^), 

(b)  ua(  • ,  •  ;x)  =  0,  uaP(  • ,  *  jx)  =  0  on  B*(-oo,oo). 

Observe  that  requirements  (a),(b),  because  of  (3.24)  and 
Theorem  1.4,  for  fixed  x€R,  uniquely  characterize  the  states 
and  f°P(  • , • ;x)  as  the  respective  solutions  of  two 
first  boundary-value  problems  for  the  region  R.  Consequently, 
the  Green’s  states  |a  and  are  also  uniquely  determined  by 
the  foregoing  definition.  With  reference  to  this  definition  we 
state 

Theorem  3.1  (Integral  representation  of  the  solution  to  the 
first  boundary-value  problem) .  Suppose 

1=  [u,e,o]€Y(H,G1,02,P).  (3.25) 

Let  /°(£,t;x)  and  |°^(S.,t;x)  be  the  Green's  states  of  the  first 
kind  for  the  region  R  and  relaxation  functions  G-^Gg.  Then»  for 
each  (x,t)  Q  R*(-oo  ,oo  ), 

ua(x,t)  *  I  [F*dua](£,t;x)dVV  -  I  [§a#du](£,t ;x)dAr  ,  (S^)1-^ 
=  J*  [F*duaP](^,t jx)dV^  -  I [SaP#du](^.,t;x)dA^.  (3.27) 

^  Here  and  in  the  sequel  we  conveniently  write  (£#d^](^,t ;x)  in 
place  of  [£(•, -jx)  *<%(•>•  ;x)  ](i>t)  if  £(  •,  •  ;x),£(  •,  •  ;x)  are 
suitable  functions  of  position  and  time .  Cf .  the  notations 
adopted  in  (l .9), (l .15; • 
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Proof .  Choose  (x,t ) £  R*(-oo  ,oo  )  and  hold  (x,t)  fixed  throughout 
the  following  argument.  Let  pQ  >  0  be  a  number  such  that 

(x)CR.  For  each  p€(0,p  ]  denote  by  Rft  the  regular  region 

~  O  °  r 

R“Qp(x)  with  the  boundary  BUZp(x).  Bearing  in  mind  Defini¬ 
tion  3.1,  as  well  as  Theorems  2.2,  2.3,  observe  that 

jV.-sx),  JoP(-.-jx)6-V{l?p.01.^).  (3.28) 

Next,  apply  the  reciprocal  theorem  (Theorem  1.5)  to  the 
pair  of  states  J,  Ja(-,-;x)  and  to  the  pair  /,  JaP(-,*jx)  in  their 
common  domain  of  regularity  Rpx(-oo  ,oo  ) .  Because  of  (b)  in 
Definition  3.1,  thi3  yields 

J  lS*duaJ(£,t;x)dAp  +  J  [F#duaKl,t;x)dVc  = 

Sp(x)  Rp 

J  [Sa*du](^,t;x)dA.  +  |  [S°#du](£,t  jxJdAj.  ,  (3.29) 

B  Sp(x) 

[S#duaP)(£,t;x)dA,  +  J  (F#du«P](i,t;x)dVr  = 

£p(x)  Rp 

J  (SaP*du ] ( J>, t j x) dA^  +  J  [Sa^du](I,t;x)dA^  .  (3-30) 

B  "  Sp(x) 

Now  proceed  to  the  limit  as  p  ->  0  in  (3.29)  and  (3.30), 
using  Lemma  3.1  and  Lemma  3.2,  respectively.1^  This  confirms 
(3. 26), (3. 27)  and  completes  the  proof. 


1  Recall  from  (3.24)  that  the  singular  part  of  Ja(*>*jx)  and  of 
i°P(*,*,x)  are  a  Kelvin-state  and  a  linear  combination  of 


doublet -states,  respectively.  On  the  other  hand,  the  regular 
parts  of  the  Green's  states  under  consideration  trivially  fail 
to  contribute  to  the  limits,  as  p  -»  0,  of  the  surface  integrals 
over  Ep(x)  in  (3. 29), (3. 30). 
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The  relevance  of  Theorem  3.1  to  the  first  boundary-value 
problem  stems  from  the  fact  that  the  right-hand  members  of  (3.26), 
(3.27)  —  apart  from  elements  of  the  Green's  states  ja ,  — 

involve  only  the  body  forces  and  surface  displacements  of  the 
state  /  .  Integral  representations  for  the  stresses  belonging 
to  1  ,  are  immediately  obtainable  from  (3.27).  If  In 
particular,  coincides  with  the  entire  space  E,  then  ?a  *  on 
Ex( -ao ,00 )xE,  according  to  Definition  3*1>  and  (3.26)  reduces  to 

uQ(x,t)  =  -I"  lj>dua)(£,t  ;x)dV£  .  (3-3l) 

E  — 

We  now  further  examine  the  Green's  states  entering 
Theorem  3.1  and  establish 

Theorem  3.2  (Symmetry  of  the  Green's  states  of  the  first  kind). 

Let  Ja(5_,tjx)  and  ]a^(§u.»t;x)  be  the  Green ' s  states  of  the  first 
kind  for  a  region  R.  Then,  for  each  (^, t,x) €  R*( -00 , 00 )xR  with 

L  /  X’ 

updL’tJ*)  *  *  (3.32) 

(3.33) 

Proof .  Choose  and  fix  (£.,t,x)  €.Rx(-co,co  )xR  with  £  ^  x.  Let 

p  >  0  be  a  number  such  that  Z  (£)  C  R,  Z_  (x)CR,  vjhile 
0  po  Fo 

Z  (OOZ  (x)  is  empty.  Then,  for  each  p£(0,p  1,  the  region 

P.  0*  “  O 

o  r  o 

Rp  -  R-Dp(D  -  5p  (x)  is  regular  and,  by  hypothesis  and 
Definition  3-1, 


562(25)/21 


41 


Je(-.-iS),  jP(-,-;i)6-Y(Kp,01,0a),  (3. 34) 

|aP<  •,  •  •.  -i)  €  •V(Ip,01,02) .  (3.35) 

To  complete  the  argument,  apply  Theorem  1.5  to  each  of  the 
pairs  of  states  appearing  in  (3.34)  and  (3.35),  respectively, 
bear  in  mind  (b)  in  Definition  3*1,  and  proceed  to  the  limit  as 
P  — >  0,  making  use  once  again  of  Lenina  3.1  and  Lemma  3.2. 

We  turn  now  to  the  second  boundary- value  problem  (surface 
tractions  prescribed  over  the  entire  boundary  for  all  time)  and 
adopt 

Definition  3.2  (Green's  states  of  the  second  kind).  We  call 
i  a  (£,  t  ;x)  =  [ua  (£,  t ;  x)  ,  £a  (£,  t ; x )  ,3°  (£,  t  ;x)  ] , 
faP(£*t;x)  “  (ia^(£»t;x),£aP(^,t;x),£aP(i,t;x)  ], 


the  Green 1 a  states  of  the  second  kind  for  a  (regular)  region  R, 
relaxation  functions  Q^Gg,  and  —  in  case  R  is  bounded  —  for 
a  (fixed  point )  x°€  R  if  and  only  if:  for  all 
(£,t,x) € "Rx(-co  ,oo  )xR  with  £  /  x#  £  /  x°. 


fa(£,tjx)  »  /a(£,tjx)  +Ja(£,tjx) 

+  c[-fa(£,t;x°)  +  ^(3.36) 

Jap(£,t;x)  =  -  |UaP(£,t;x)  +/Pa(£,t;x)]  +  JaP(£,t;x), 

where  Ja(£,tjx),Ja^(£,t;x),  and  7°  (£***£)  are  the  normalized 
Kelvin-state,  the  doublet-state,  and  the  state  appropriate  to 
a  center  of  rotation,  all  corresponding  to  G1,Gg,  while  c=l 
when  R  _ls  bounded  and  c=0  when  R  is  unbounded :  further. 


562(25)/21 


42 


j°(I,t;x)  =  [u°(^tjx),ia(^,tjx),£a(i,tjx)], 

Ja^(Ji»tjx)  -  fuaP(^Jt;x)»Ia^(Ijtjx)>3aP(^>t]x)  ], 
for  each  x£R>  are  states  such  that 

(a)  J°( ' >  *  ;x)*  Ja^(  *,  -jx)  €  *V(T?, O^G^), 

(b)  Sa(*,-jx)  =  0,  SaP(*,*;x)  =  0  on  Bx( -00,00), 
and 

(c)  when  R  Is  bounded, 

ua(x°,  •  ;x)  =wa(x°J-;x)  =0  on  (-00,00), 
uaP(x°,*;x)  *  (x°,  •  ;x)  =  0  on  (-00,00), 

In  which  ilia ,  ZP& are  the  respective  rotation  vectors  of  Ja,  J°P. 

Definition  3.2  is,  to  an  extent,  analogous  to  Defini¬ 
tion  3.1,  The  regular  parts  of  the  Green's  states  in 
Definition  3.1  are  uniquely  characterized  as  solutions  to  first 
boundary-value  problems  for  the  region  R.  In  contrast,  Ja(*,*jx) 
and  JaP{ • , • jx) ,  for  fixed  x  €  R,  are  in  the  present  instance  each 
the  solution  of  a  second  boundary-value  problem  for  R,  as  is 
apparent  from  (a),  (b),  and  (3*36).  According  to  Wieorem  1.4, 
the  solution  to  such  a  problem  is  unique  —  except  for  an 
arbitrary  additive  (infinitesimal)  rigid  motion  of  the  entire 
body,  when  R  is  a  bounded  region. ^  Condition  (c)  in  Defini¬ 
tion  3.2  serves  to  eliminate  this  indeterminacy  from  uP  and  tPP . 

if  r  is  unbounded,  the  additive  rigid  motion  is  precluded  by 
the  regularity  condition  (c)  in  Definition  1.3,  which  requires 
the  displacements  to  vanish  at  infinity.  Note  also  that  in 
the  flrBt  boundary- value  problem  the  additive  rigid  motion  is 
precluded  by  the  boundary  conditions,  regardless  of  whether 
or  not  R  is  bounded. 
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Thus,  for  a  given  x°*  the  states  ja,JaP  are  fully  unique  and 
hence,  because  of  (3.36),  the  same  is  true  of  the  Green's  states 
Ja,ia^  of  the  second  kind.  Further,  when  R  is  unbounded  (c=0), 
both  fa  and  jaP  are  entirely  independent  of  x°. 

If  R  is  bounded  (c=l),  the  singular  part  of  fa  in  (3-36) 
is  considerably  more  complicated  than  its  counterpart  in  (3.24) 
of  Definition  3.1  for  reasons  to  be  made  clear  presently. 

Indeed,  a  necessary  condition  that  the  second  boundary-value 
problem  characterizing  fa  possess  a  solution  when  R  is  finite, 
is  that  for  each  x£R, 

J  ±>a(iL,  *;x)dA£  =  0,  J  £.ASa(£,  *  jx)dAj.  =  0  on  (-00,00 ), 

B  B  (3.37) 

i.e.  that  the  surface  tractions  governing  Ja(*,*;x)  be  self- 
equilibrated The  requirement  (3-37),  in  turn,  because  of  (b) 
and  (3.36),  implies  that  the  system  of  singularities  involved 
in  the  singular  part  of  fa  must  be  self -equilibrated .  The 
supplementary  singular  part  of  fa,  which  carries  the  multiplier 
c  and  whose  singularities  are  located  at  K  -  x°*  serves  the 
purpose  of  assuring  the  self-equilibrance  of  the  complete  system 
of  singularities  at  £  =  x  and  K  =  x° .  This  claim  is  readily 
verified  with  the  aid  of  (b),(c)  in  Theorem  2.2  and  (2.42). 

Note  from  (a)  that  must  meet  the  equilibrium 

equations  in  the  absence  of  body  forces . 
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Conditions  (3.37)  are  no  longer  necessary  for  the  exist¬ 
ence  of  Ja  when  R  is  unbounded,  in  which  case  c=0  and  the 
singular  part  of  is  merely  a  normalized  Kelvin-state  corre¬ 
sponding  to  a  concentrated  load  at  £  =  x.  Finally,  we  observe 
that  the  pair  of  singularities  at  £  =  x  in  the  singular  part  of 

AdB 

j  is  always  self-equilibrated,  regardless  of  whether  or  not  R 
is  finite,  as  is  immediate  from  the  second  of  (3.36)  and  (b),(c) 
in  Theorem  2.3*  We  may  now  proceed  to 

Theorem  3.3  (Integral  representation  of  the  solution  to  the 
second  boundary-value  problem) .  Suppose 

d]  €7(^,0^ G2,F)  (3.38) 

and,  if  R  is  bounded, 

u(x°,  *)  51  “(x°»  •)  =  0  on  (-co, oo),  (3.39) 

where  is.  the  rotation  vector  belonging  to  J  ,  while  x°G  R. 

Ja(£.,tjx)  and  JaP(£.,t;x)  be  the  Green 's  states  of  the  second 
Wind  for  the  region  R,  the  relaxation  functions  0^,0^  and  —  In 
case  R  Is  bounded  —  for  x° .  Then,  for  each  (x.t)  £  Rxf-m  ) . 

V^t)  =  J  [F*djua ] (£_, t ;x)dVp  +  J[S#dua](£,t;x)dAr  ,  (3.40) 

R  B  jl 

=  ^  [P*dMa^Ki>t;x)dV^  +  Jf^du®^  ](£,  t  ;x)dA^  .  (3 .41) 

Proof.  Note  that  if  i  satisfies  (3-38)  with  R  bounded,  it  can 
always  be  made  to  meet  (3-39)  as  well  by  addition  to  u  of  a  rigid 
displacement  history.  The  indentlty  (3.4l)  and,  if  R  is  unbounded, 
also  (3.40),  may  be  confirmed  by  an  argument  strictly  analogous 
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Conditions  (3.37)  are  no  longer  necessary  for  the  exist¬ 
ence  of  Ja  when  R  is  unbounded,  in  which  case  c=0  and  the 
singular  part  of  |a  is  merely  a  normalized  Kelvin-state  corre¬ 
sponding  to  a  concentrated  load  at  £  =  x.  Finally,  we  observe 
that  the  pair  of  singularities  at  =  x  in  the  singular  part  of 
is  always  self-equilibrated,  regardless  of  whether  or  not  R 
is  finite,  as  is  immediate  from  the  second  of  (3.3*5)  and  (b),(c) 
in  lheorem  2.3*  We  may  now  proceed  to 

Theorem  3*3  (Integral  representation  of  the  solution  to  the 
second  boundary-value  problem) .  Suppose 

[u  ,e,d]  £  YO^Q^O^F)  (3.38) 

and,  if  R  is  bounded, 

u(x°,  • )  *  w(x°,  *)  “  0  on  (-co, oo),  (3-39) 

where  ^  is  the  rotation  vector  belonging  to  J  ,  while  x°€  R. 

Let  ja(£,tjx)  and  jaP(£,t;x)  be  the  Green 's  states  of  the  second 
kind  for  the  region  R,  the  relaxation  functions  G-j^Qg*  and  —  in 
case  R  is  bounded  —  for  x° .  Then,  for  each  ( x, t ) £  Rx( -oo ,  oo  ) , 

VLjfxjt)  =  J  [F#djQa](l,t;x)dVr  +  J[S*dua](£, t ;x)dAj.  ,  (3.40) 

R  *  B  i 

eap(x,t)  =  J  (F*duaP](^,t;x)dV^  +  JfSedjj®? ](£, t ;x)dA^ .  (3.41) 

Proof.  Note  that  if  S  satisfies  (3.38)  with  R  bounded,  it  can 
always  be  made  to  meet  (3-39)  as  well  by  addition  to  u  of  a  rigid 
displacement  history.  The  indentity  (3.41)  and,  if  R  is  unbounded, 
also  (3.40),  may  be  confirmed  by  an  argument  strictly  analogous 
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to  that  employed  in  the  proof  Theorem  3.1 •  We  therefore  give  a 
detailed  derivation  merely  for  (3.40),  assuming  R  to  be  bounded. 

Consider  first  the  special  case  in  which  x  =  x°  and  con¬ 
clude  from  (3*36)  that  then 

Ja(-i*;x0)  =Ja(*>,;x°)  on  T?x(-oo  ,00 ) .  (3-42) 

But  (3.42),  together  with  (a),(b),(c)  of  Definition  3.2  and 
Theorem  1.4,  implies  that  J^*,*;*0)  is  the  null-state.  Hence, 
and  in  view  of  the  first  of  (3*39)*  the  identity  (3-40)  is 
trivially  met  when  x  =  x°* 

Next,  hold  (x,t)  €  Rx(-oo  ,00 )  fixed  and  assume  x  /  x°. 

Let  p  >  0  be  a  number  such  that  2  (x)CR,  2  (x°)CR,  while 

0  po  po 

2  fe)n£0  (x°)  is  empty.  Then,  for  each  p€(0,p  ],  the  region 
po  po 

Rp  ■  R-C^,(x)-Qp(x°)  with  the  boundary  BUEp(x)u  2p(x°)  is  regular 
and,  by  hypothesis  and  Definition  3.2, 

?a(*,-;x)fe'V(H,,o1,o2).  (3.43) 

Now  apply  the  reciprocal  theorem  (Theorem  1.5)  to  the  pair  of 
states  i  and  ( • , • jx)  in  their  oommon  domain  of  regularity 
lLx(-oo,oo).  Because  of  (b)  in  Definition  3.2  one  thus  arrives  at 


J  [P#du°  ](5.,t;x)dVr  + 

J  [S*du°](^,t;x)dAr  + 

B  i 

J  [  ]  ( i,  t ;  x )  dkr  + 

J  (S#dua](^,t;x)dA^  = 

2p(x) 

2p(x°) 

[§°*du](£.,t jx)dAr  + 

J  [Sa*du)(£,t;x)dAr  . 

y*> 

2p(x°) 

(3.44) 


562(25)/21 


46 


Proceed  to  the  limit  as  p  — *  0  in  (3.44),  taking  into  account  the 
nature  and  location  of  the  singularities  entering  the  right-hand 
member  in  the  first  of  (3*36)  and  using  Lemma  3.1*  Corollary  3.1, 
as  well  as  (3.39).  This  yields  the  desired  result  (3.40). 

It  is  instructive  to  examine  the  influence  of  the  choice 
of  the  fixed  point  x°  in  Theorem  3-3.  A  change  in  this  choice 
evidently  affects  u  merely  within  an  additive  rigid  displacement 
history  and  therefore  leaves  £  unaltered.  On  the  other  hand, 
such  a  change  alters  the  basic  structure  of  Qa  (whose  singular 
part  depends  upon  the  location  of  x°)  but  results  only  in  the 
addition  of  a  rigid  displacement  to  uaP,  which  is  easily  seen  to 
have  no  effect  on  the  right-hand  member  of  (3.4l). 

The  Green's  state  faP  of  the  second  kind  conforms  to  the 
symmetry  relation  (3.33)  In  Theorem  3.2,  as  may  be  verified  by 
precisely  the  same  scheme  used  to  prove  (3-33)  originally.  In 
contrast,  Ja  of  Definition  3.2  Is  found  not  to  obey  (3.32)  because 
of  the  asymmetric  manner  in  which  £  and  x  enter  the  supplementary 
singular  part  of  ja  in  the  first  of  (3.36). 

Iheorem  3*3  owes  its  importance,  as  far  as  the  second 
boundary-value  problem  is  concerned,  to  the  fact  that  the  right- 
hand  members  of  (3.40), (3.41)  —  apart  from  elements  of  the 
Green's  states  of  the  second  kind  —  involve  exclusively  the  body 
forces  and  surface  tractions  appropriate  to  J .  Analogous  Integral 
representations  for  the  stresses  belonging  to  /  follow  at  once 
from  (3.41)  with  the  aid  of  (l .19), (1 .20) . 
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We  have  so  far  confined  our  attention  to  integral  repre¬ 
sentations  appropriate  to  the  first  and  second  fundamental 
boundary-value  problems.  An  integral -representation  theorem  for 
the  mixed  problem,  which  contains  Theorem  3.1  and  Theorem  3-3  as 
special  cases,  may  be  deduced  by  obvious  similar  means  and  with 
but  superficial  complications. 

Theorem  3*3  may  be  used  as  a  basis  for  a  mathematically 
precise  and  physically  meaningful  definition  of  the  notion  of  a 
concentrated  surface  load .  The  latter  concept  may  be  defined 
through  a  limit  process  applied  to  a  sequence  of  regular  visco- 

17 

elastic  states  that  corresponds  to  distributed  surface  tractions, 
in  analogy  to  the  limit  treatment  of  internal  concentrated  loads 
contained  in  Section  2.  We  shall,  however,  not  pursue  this  issue 
further  and  shall  turn  instead  to  a  more  basic  application  of  the 
integral  representation  for  the  solution  to  the  second  boundary- 
value  problem  supplied  by  Theorem  3.3. 


17 


See  [8],  Section  7,  for  a  detailed  analysis  of  concentrated 
surface  loads  in  elastostatics . 
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4.  Saint-Venant 's  principle  for  viscoelastic  solids. 

Saint -Venant 's  principle  in  the  classical  equilibrium 
theory  of  elastic  solids  was  originally  introduced  by 
Saint-Venant  [153  in  connection  with  —  and  with  limitation  to  — 
the  problem  of  extension,  torsion,  and  flexure  of  prismatic  or 
cylindrical  bodies .  The  earliest  universal  statement  of  the 
principle  is  apparently  due  to  Boussinesq  [l6](p.298),  whose 
formulation  has  since  become  traditional.  Love  [5](p*132), 
adopting  Boussinesq' s  version,  states  the  principle  as  follows: 
"...  the  strains  that  are  produced  in  a  body  by  the  application,  to 
a  small  part  of  its  surface,  of  a  system  of  forces  statically 
equivalent  to  zero  force  and  zero  couple,  are  of  negligible 
magnitude  at  distances  which  are  large  compared  with  the  linear 
dimensions  of  the  part."  Hie  far-reaching  importance  of  the 
principle  stems,  of  course,  from  the  fact  that  it  presumably 
entitles  one  to  relax  the  boundary  conditions  in  the  second 
boundary-value  problem  of  elastostatlcs  and  to  replace  the  given 
surface  tractions,  at  least  in  part,  by  statically  equivalent— 
but  analytically  more  manageable  —  loadings. 

As  was  first  pointed  out  by  von  Mises  [6],  the  conven¬ 
tional  statement  of  Saint-Venant 's  principle  is  in  need  of 
clarification.  Von  Mises  observed  that  the  sentence  cited 
involves  a  dual  comparison,  which  is  not  fully  made  explicit. 

Thus,  the  statement  in  question  asserts  that  the  strains  due  to 
self -equilibrated  loads  are  "negligible"  at  distances  which  are 
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large  compared  to  the  size  of  the  load-region,  the  tacit  implica¬ 
tion  being  that  these  strains  are  small  compared  to  strains  pro¬ 
duced  by  loads  which  are  not  statically  equivalent  to  null . 

Von  Mises  noted  further  that  one  cannot  in  general  meaningfully 
speak  of  strains  "produced"  by  non-equilibrated  loads  since  the 
equilibrium  problem  for  given  surface  tractions  and  for  a  bounded 
region  (in  the  absence  of  body  forces)  has  no  solution  unless  the 
tractions  acting  on  the  entire  boundary  conform  to  equilibrium. 

He  concluded  that  for  a  finite  body  it  is  necessary  to  consider 
the  Joint  effect  of  the  tractions  applied  to  several  (at  least 
two)  distinct  portions  of  the  boundary. 

Two  additional  criticisms  of  the  traditional  version  of 
Saint - Venant 's  principle  are  equally  self-evident.  Trivially, 
the  strains  arising  from  loads  applied  to  a  finite  part  of  the 
surface  of  an  unbounded  elastic  body  are  arbitrarily  small  at 
points  sufficiently  far  removed  from  the  region  of  load  applica¬ 
tion,  regardless  of  whether  or  not  the  loading  is  equilibrated. 

On  the  other  hand,  it  follows  from  the  superposition  principle 
of  linear  elasticity  theory  that  the  strains  and  stresses  at  a 
fixed  point  of  an  elastic  body,  produced  by  equilibrated  surface 
loads,  may  be  made  as  large  as  one  pleases  by  choosing  the 
magnitude  of  the  loads  sufficiently  large,  regardless  of  the 
size  of  the  region  of  load  application. 

The  foregoing  observations  suggest  an  interpretation  of 
the  conventional  statement  of  Saint- Venant 's  principle  that  may 


562(25)/21 


50 


lR 

roughly  be  phrased  as  follows  :  Let  the  loads  acting  on  an 
elastic  body  be  confined  to  several  distinct  portions  of  its 
boundary,  each  lying  within  a  sphere  of  radius  p,  and  suppose  the 
loads  remain  bounded  as  p  — >  0.  Then  the  strains  at  a  fixed 
interior  point  of  the  body  are  of  a  smaller  order  of  magnitude 
in  p ,  as  p  — >  0,  when  the  tractions  on  each  load  region  are  self- 
equilibrated  than  when  they  are  not. 

That  this  is  the  meaning  intended  by  Boussinesq  [16]  is 
apparent  from  his  own  efforts  to  justify  the  principle.  With 
this  objective  in  mind  he  examined  the  strains  at  an  interior 
point  of  a  semi-infinite  elastic  body  that  is  subjected  to  a  set 
of  concentrated  loads  applied  normal  to  its  plane  boundary. 
Assuming  the  points  of  application  of  the  loads  to  lie  within  a 
sphere  of  radius  p,  he  then  showed  that  the  order  of  magnitude 
of  the  strains  in  question  is  p  if  the  resultant  force  is  zero, 

p 

and  p  when  the  resultant  moment  also  vanishes. 

Von  Mises  [6]  demonstrated  with  the  aid  of  two  examples, 
which  involve  tangential  as  well  as  normal  surface  loads,  that 
the  traditional  version  of  Saint -Venant 's  principle  (interpreted 
as  above)  requires  amendment  in  order  to  be  generally  valid. 
Guided  by  these  counter-examples  he  conjectured  a  modified 
Saint-Venant  principle  which  was  later  on  formulated  and  proved 
in  [7],  It  is  this  modified  principle  which  we  now  seek  to 
extend  to  viscoelastic  solids . 

This  interpretation  is  taken  from  [73;  it  is  an  elaboration 
of  von  Mises'  [6]  earlier  interpretation. 
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In  order  to  avoid  an  unduly  lengthy  statement  of  the 
theorem  to  be  established  we  introduce  two  preliminary  definitions 
that  will  permit  us  to  phrase  the  underlying  hypotheses  concisely. 
In  this  connection  we  recall  our  previous  agreement  to  the  effect 
that  Qp(x°)  always  stands  for  an  open  sphere  of  radius  p 
centered  at  x°. 

Definition  4.1  (A  set  of  families  of  contracting  load  regions). 

We  say  that  Ap  (0  <  p  <  pQ;  n=l,2, . . .N)  is  a  3et  of  N  families 

of  load  regions  on  the  boundary  B  of  a  (regular)  region  R,  which 

c 

contract  to  N  points  B,  if  and  only  if  for  every  p  £ (0,pQ) 

and  every  n(n=l,2, . . .N) : 

_  o 

(a)  Ap  =  Qp(jj_n)r*B  and  this  intersection  is  connected; 

(b)  0  A^  =  0  for  m  /  n  (m=l,2,  .N) ; 

(c)  A^1  C  nn  i  where  nn  (the  "embedding  region"  of  Ap) 

13  a  subset  of  B  and  the  position  vector  £_  of  nn  admits  the 
parametrlzatlon 

L  =  In(i)  for  a  =  (a-^a^eD11,  (4.1) 

Dn  being  an  ojcen,  bounded,  simply-connected  plane  region,  while 

fn  £  C2(Dn),  fn(0)  =  ln,  f^A  f^2  fi  0  on  Dn.  (4.2)19 

A  typical  member  of  the  set  A^  (0  <  p  <  pQj  n=l,2, ...N), 
for  fixed  p  and  n,  is  shown  in  the  accompanying  figure,  in  which 

15 - - - T“ 

We  use  the  notation  f "  p(a)  =  35-  fn(a)  (p=l,2) .  No  summation 
with  respect  to  n  is  intended  inP(4.2). 
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Q  and  Q  designate  the  endpoints  of  the  position  vectors  and 

K,  respectively,  whereas  x  is  the  position  vector  of  a  point  P 

in  R.  Note  that  according  to  the  first  of  (4.2)  each  subregion 

nn  of  B  possesses  continuous  curvatures.  Therefore  the  points 
°n 

V  (n=l,2, . . .N)  are  necessarily  distinct  regular  points  of  B; 

moreover,  the  boundary  B  —  within  each  of  the  N  embedding 

regions  —  is,  by  implication,  required  to  exhibit  a  higher 

degree  of  smoothness  than  that  automatically  assured  by  the 

assumption  that  R  is  a  regular  region  of  space.  The  second  of 

o 

(4.2)  asserts  merely  that  Kf1  is  the  image  under  fn  of  the  origin 

of  the  parameter-plane  ( "a -plane " ) .  Finally,  the  third  of  (4.2) 

is  equivalent  to  the  condition  that  at  each  point  of  Dn  at  least 

one  of  the  Jacoblans  of  the  mapping  £n  fails  to  vanish,  so  that 

fn  defines  a  regular  curvilinear  coordinate -net  on  IIn. 

Definition  4.2  (Associated  family  of  viscoelastic  states).  Let 

Ap  (0  <  p  <  pQ ;  n=l,2, . . .N)  be  a  set  of  N  families  of  load  regions 

on  the  boundary  B  of  a  region  R,  which  contract  to  N  points 
o 

j»n£B.  We  say  that 

<f(x*t,p)  =  [u(x, t,p),.e(x,t,p),jj(x,t,p)  ]  (0  <  p  <  pQ)(4.3) 

is  a  family  of  viscoelastic  states  on  Ex( -00 , 00 )  corresponding  to 
loads  on  a£  if  and  only  if  for  each  p  € (0,pQ) : 

(a)  J  ( •,  *,p)  €'V(l,a1,a2) ; 

N 

(b)  S(-,*,p)  =  0  on  (B  -  U  A*:)x(-co  ,co ) ; 

n=l  p 

(c)  when  R  _is  bounded. 
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Q  and  Q  designate  the  endpoints  of  the  position  vectors  Jj”  and 

K,  respectively,  whereas  x  is  the  position  vector  of  a  point  P 

in  R.  Note  that  according  to  the  first  of  (4.2)  each  subregion 

nn  of  B  possesses  continuous  curvatures.  Therefore  the  points 
°n 

V  (n=l,2, . . .N)  are  necessarily  distinct  regular  points  of  B; 

moreover,  the  boundary  B  —  within  each  of  the  N  embedding 

regions  —  is,  by  implication,  required  to  exhibit  a  higher 

degree  of  smoothness  than  that  automatically  assured  by  the 

assumption  that  R  is  a  regular  region  of  space.  The  second  of 

o  « 

(4.2)  asserts  merely  that  V1  is  the  image  under  f“  of  the  origin 

of  the  parameter-plane  ("a-plane").  Finally,  the  third  of  (4.2) 

is  equivalent  to  the  condition  that  at  each  point  of  Dn  at  least 

one  of  the  Jacobians  of  the  mapping  fn  fails  to  vanish,  so  that 

fn  defines  a  regular  curvilinear  coordinate-net  on  IIn. 

Definition  4.2  (Associated  family  of  viscoelastic  states).  Let 

Ap  (0  <  p  <  pQ;  n=l,2, . . .N)  be  a  set  of  N  families  of  load  regions 

on  the  boundary  B  of  a  region  R,  which  contract  to  N  points 
o 

B.  We  Bay  that 

^(x,t,p)  =  [u(x,t,p),£(x,t,p),o(x,t,p) )  (0  <  p  <  pq)(4.3) 

is  a  family  of  viscoelastic  states  on  Kx( -00 , 00 )  corresponding  to 
loads  on  if  and  only  if  for  each  p  € (0,pQ) : 

(a)  )£'V(E,Q1,G2); 

N 

(b)  S(-,-,p)  =  0  on  (B  -  U  A^)x(-oo,co ); 

n=l  p 

(c)  when  R  lj»  bounded. 
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u(x°, *,p)  =  w(x°>  '•> p)  =  0  on  (-00, oo), 
where  w  is  the  rotation  vector  of  4  and  x°€  R; 


(d)  |Sj  Is  uniformly  bounded  on  Bx(-oo,t ]x(0, pQ)  for 
every  t  € [0,oo) . 

Furthermore  we  adopt  the  notation 


Ln(t,p) 


tf(t,p) 


—  J  S(£,t,p)dA  , 

=  J  lAS(^,t,p)dA  , 


>  (4.4) 


whence  Ln(t,p)  and  ff*(t,p)  stand  for  the  resultant  force  and  the 
resultant  moment  about  the  origin,  of  the  tractions  S( *,t,p) 
acting  on  a”  . 

Assumption  (b)  requires  the  surface  tractions  S(*,t,p) 
of  the  state  /  (*,t,p)  to  vanish  on  the  complement  with  respect  to 
B  of  the  union  of  the  load  regions  (n»l,2, . . .N),  for  all  time 
and  every  p  € (0,pQ) .  It  Is  worth  noting  that  no  regularity 
restrictions  other  than  (d)  are  placed  on  the  family  of  states 
4  as  far  as  Its  dependence  upon  the  parameter  p  Is  concerned. 
Indeed,  but  for  notatlonal  complications.  It  would  have  been 
equally  adequate  for  our  purposes  to  Introduce  a  set  of  sequences 
of  contracting  load  regions  and  an  associated  sequence  of  visco¬ 
elastic  states  corresponding  to  loads  on  these  subregions  of  the 
boundary.  We  may  now  proceed  to  formulate 
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Theorem  4.1  (A  Saint -Venant  principle  for  viscoelastic  solids). 

Let  Ap  (0  <  p  <  pQj  n=l,2, . . .N)  be  a  set  of  N  families  of  load 

regions  on  the  boundary  B  of  a  (regular)  region  R,  which  contract 

o  ~ 

to  N  points  £.  €  B,  and  assume  N  >  2  if  R  is  bounded.  Let 
/ (0  <  p  <  pQ)  be  a  family  of  viscoelastic  states  on 
ftx(-oojoo)  corresponding  to  loads  on  Ap  .  Let  xfcR  and  t€ (-00,00) 
Further,  assume  the  existence  of  the  Qreen 1 s  states  of  the  second 
kind  for  the  region  R  (Definition  3.2) .  Then,  uniformly  on  (-oo,t  ] 


n(Xi  •  >P )  =  0(p6),  e(x, -,p)  =  0(p5),  d(x, -,p)  *  0(p6)  as  p  0, 

(4.5) 

where  6  =2 .  Moreover : 


(a)  6  *  3  If 

Ln  *  0  on  (-oo,t]x(0,po),  (n*l,2, . . .N); 

(b)  6  *  4  if 

Ln  -  0  ,  J  *,-) dA  -  0 

An 

AP 

on  (-oo,t)x(0,po),  (n»l,2, . ..N); 


(4.6) 


?  (4.7) 


(c)  6  «  4  if 

Ln  =  0,  (f  «  0  on  (-oo,t]x(0,pQ),  (n=l,2, . . .N),  (4.8) 

provided 

S(L,*c,p)  =  ®n(L*T*P)i^1(T)  (no  sum)  (4.9) 

for  all  (L,*c,p)  €  ApX(-oo,t]x(o,po)  and  n-1,2,... N.  Here  4>n  is 
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scalar-valued .  while  kf1  Is  continuous  with  kn(x)  a  unit  vector 
such  that 

kn-?n  ^  0  on  (-00, t]  (n=l,2, . .  .N),  (4.10) 

On  °n 

V  being  the  unit  normal  of  B  at  V. 


Proof.  By  hypothesis,  (a)  In  Definition  4.2,  and  Definition  1.3, 

/  (x,t,p)  Is  the  null  state  for  each  (x,t,p)€  Rx(-oo  j,0)*(0, p  ); 

”  ~  o 

accordingly  the  conclusion  is  trivial  when  -oo  <  t  <  0.  Thus 
choose  (x,t) €  Rx[0,oo )  and  hold  (x,t)  fixed  for  the  remainder  of 
the  argument . 

Prom  (a),(b),(c)  in  Definition  4.2,  (b)  in  Definition  4.1, 
Theorem  3-3,  and  (1.19),(1.20)  follows,  for  all 
(**,P)  €(-co,co)x(o,po). 


Ui(x,x,p)  ■  2  u"(x,x,p), 

n»=l 

N 

®  i  i  (i»"^  »P )  =  £  £j  *(x,t,p), 

n»=l  1J 

N 


*> 


r  (4.n) 


where,  for  each  n  (n*l,2, . . ,N), 

“  J  [S*du1](£,‘t;x,p)dAE 
AP  A 

1(x,'c,p)  =  f  [S*du1,J)(^,xjx,p)dA£  , 

A?  1 

oJj(£,x,p)  =  [eJ^C^Hx.x.p) 

+  3  6iJ^£kk*c^°2“(V^-,x,P)  * 


>  (4.12) 


20 


W 


We  write  I Swd^1  ]  (£.,  t  ;x,p  )  to  denote  the  convolution-value 

[S(  *,  *,p)*du1(  •,  *jx))(t,t),  etc.  See  also  (l .9), (1 .15), 
well  as  Footnote  Ifc.  13.  v  37 ' 


as 
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Here  u^g^xjx)  and  ul4*(£yc;x)  are  the  displacement  field  histories 
appropriate  to  the  Green's  states  of  the  second  kind  for  the 
region  R,  the  relaxation  functions  0^,Qg,  and  —  in  case  R  is 
bounded  —  for  x°  (See  Definitions  3.2,  4.2). 

In  view  of  (4.11), (4.12)  it  is  natural  to  call  uJ,eJj, 
andtJ^j  (n=l,2,...N)  the  displacement,  strain,  and  stress  contri¬ 
butions  arising  from  the  loading  on  the  n-th  family  of  load 
regions.  We  mention  parenthetically,  however,  that  when  R  is 
bounded  the  "contribution-state"  /n  =  [un,e_n,d.n]  (n=l,2,...N) 
possesses  an  independent  physical  significance  as  the  solution 
of  a  second  boundary- value  problem  if  and  only  if  the  tractions 
on  each  individual  family  of  load  regions  are  permanently  self- 
equilibrated,  i  .e.  Ln  »  Mf1  *  0  on  (-00 ,00  )x(0,pQ)  for  n*=l,2,  ...N. 
In  this  case,  or  when  R  is  unbounded,  one  has  for  each  pG(0,po) 
and  each  n  (n=l,2, . . .N), 

with 

sn( *, *,p)  =  S( •, -,p)  on  Ap  , 

Sn(.,.,p)  =  0  on  B-A”  , 

which  characterize  Jn  uniquely.  Also,  clearly,  when  R  is  bounded 
either  the  loading  on  each  family  of  load  regions  is  permanently 
self -equilibrated  or  there  exist  at  least  two  such  families  for 
which  this  is  not  true. 


(4.13) 

* 

v  (4.14) 
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We  now  examine  a  typical  displacement  contribution 
un(x,T,p)  with  a  view  toward  establishing  the  order  of  magnitude 
of  u(x,x,p)  as  p  — >  0.  For  this  purpose  hold  n  (n=l,2, . . .N) 
fixed  and  observe  on  the  basis  of  (c)  in  Definition  4.1  that, 
for  all  (5_,t  )€  Ilny( -00,00  ), 

u^L/tjx)  =  ju1(f^n(q_) ;x)  =  £ln(a  ,x  ;x) ,  o*(a1,o2),  (4.15) 

where  the  functions  &in(*,*,*x)  are  defined  on  Dnx(-oo,co)  and 
evidently 

fiin( * , * ,*x)  =  0  on  Dnx(-oo,0).  (4.16)21 

As  is  clear  from  the  behavior  of  u1(*,,;x)  on  Bx[0,co)  implied 
by  Definition  3-2  and  the  smoothness  of  fn  stipulated  in  the 
first  of  (4.2),  the  functions 

£ln(o.,T  ;x)  s  £in(q_,x;x)  (4.17) 

exist  and  are  continuous  for  all  (a,x)  €Dnx[0,oo ) ;  furthermore, 
throughout  this  domain  £in(a  ,x  ;x)  and  £in(a ,t ;x)  are  twice 
continuously  differentiable  with  respect  to  0i>a2- 

In  view  of  the  preceding  observations,  £in(*,x;x)  for 
each  x  €  ( -oo  ,  co  )  and  all  q^D11  admits  the  Taylor  expansion 

£ln(atx;x)  =  £in(x  jx)  +  |*£(x  ;x)ap  +  0ln(aL,x  ;x)  (4.18)22 

where  we  have  used  the  notation 

- - 

Recall  that  x€R  is  fixed. 

pp 

Here  and  in  the  sequel  summation  with  respect  to  p  (p=l,2) 
is  implied  when  p  is  a  repeated  index. 
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£ln(x;x)  =  £ln(0,Tjx), 

fi|p(Ti£)  5  3i- filn(^T;2)|(o1=a2=0)  (p*1-2)- 


(4.19) 


The  remainder  0ln(a,T;x)  in  (4.18)  evidently  possesses  the  same 
degree  of  smoothness  for  all  (ct,T) £  Dnx[0,oo )  as  does  £  (s.i'tjx). 
Furthermore, 

Qln(*,*;x)  =  0  on  dM-co^O),  (4.20) 

whereas 

0ln(a ,  •  ;x)  =  0(a2),  0la(a, • jx)  =  0(a2)  as  a  =  |a|  -»  0, 

(4.21) 

uniformly  on  (-oo,t]. 

Before  continuing  the  argument  we  note  from  (4.l),(4.2) 
that  the  mapping  fn  is  one-to-one  in  a  neighborhood  of  the  origin 
of  the  parameter-plane.  Indeed,  (4.l),(4.2)  together  with  (a) 
in  Definition  4.1  imply  that  there  exists  a  number  p-^0  <  px  <  pQ) 
and  a  function  £n  mapping  Ap  onto  a  neighborhood  of  the  origin 

of  the  a-plane,  i.e. 


£  =  £n(I), 


(4.22) 


where  £n  is  independent  of  one  of  the  components  of  £  and  is 
continuously  differentiable  with  respect  to  the  remaining  two. 

In  order  to  avoid  cumbersome  notation  and  since  n  is  being  held 
fast,  we  shall  hereafter  write  a  (£)  in  place  of  £.n(£.)  •  follows 
from  the  regularity  of  the  inverse  mapping  under  consideration  that 
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max  |a(£jl  =  0(p), 

Ap 


0(P2)  as  p  0. 


(4.23) 


Now  substitute  from  (4.15)  Into  the  first  of  (4.12)  and 
use  (4.18), (4.22) .  After  a  brief  computation  involving  permissi¬ 
ble  reversals  in  the  order  of  the  processes  of  convolution  and 
surface  integration  one  thus  obtains,  for  all  (^,p  )  €  (-oo,t  ]x(0,p^). 


uj(x,x,p)  =  I*n(x,x,p)  +  l|n(x,T>P)  +  I*n(x,x,p),  (4.24) 


-in 


in/ 


where 


xin(i*TiP)  =  lj  S(L>  -,P  )dA#d£in(  •  ;x)  ](t)  , 

Apn 

=  (J  £(£.,  ‘,p)o _(L)dA*d£!^(  *;x)  ](t)  ,  l  (4.25) 
Apn  P  'P 

lin(x,T,p)  =|  [S(£,  •,p)«d9ln(a(i),-;x)J(T)dAj  . 


,°in/ 


Also,  by  (4.16), (4.19)  and  (4.20), (4.25), 


4n(x*’>’)=0  on  (-oo,0)v(0,p1).  (4.26) 

Our  next  task  consists  in  estimating  the  order  of 
magnitude  of  I^n(x,T,p)  as  p  — >  0.  Consider  first  I^n  in  (4.25) 
and  observe  with  the  aid  of  (I.15)  and  (3-5)  that,  for  all 
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I  )  t  < 


[  max  |S(-, *,p ) J  ][  max  | ©in(a(£) ,0jx) | 
ApX[0,t]  i^Ap 


+  t  max  |eln(a(L)i^;x)N  J  dA  .  (4.27) 

(C/c)  €  A£x[o,t]  Ap 

By  (4.27),  (d)  in  Definition  4.2,  (4.23),  and  (4.21),  there  exist 
constants  p2(0  <  P2  <  P^)  and  C  such  that,  for  all 
(^,P)€  [0,t)x(0,p2), 

Il3n(x^,p)|  <  Cp4  +  Ctp4.  (4.28) 


This  conclusion,  together  with  (4.26),  assures  that 

l3n(x>‘*P)  =  0(p4)  as  p  — >  0,  (4.29) 


uniformly  on  (-oo,tJ.  Proceeding  similarly  with  the  first  two 
of  (4.25)  —  bearing  in  mind  (4.17), (4.19),  and  the  continuity  of 
£ln(0 , ’jx)  £jp(°i'jx)  on  [0,t]  —  one  arrives  at 


^n(x*-*P)  =  0(p2)  as  p  -»  0  , 
lf(x,  ',p)  =  0(p3)  as  p  -»  0  , 


j-  (4.30) 


uniformly  on  (-oo,tj.  But  (4. 29) #(4.30),  because  of  (4.24)  and 
the  first  of  (4.11),  imply  the  first  of  (4.5)  with  &=2. 

We  now  turn  to  the  remaining  assertions  concerning  the 
order  of  magnitude  of  the  displacements  as  p  — >  0,  which  pre¬ 
suppose  various  restrictions  upon  the  loading  beyond  those 
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implied  by  Definition  4.2.  Thus  suppose  (4.6)  holds  so  that  the 
resultant  force  of  the  loading  on  each  family  of  load  regions 
vanishes  up  to  the  Instant  t.  In  this  case,  I^n(x, *,•)  vanishes 
on  (-00  ,t ]x(o,p1)  according  to  (4.4), (4.25),  and  hence  (4.24), 
(4.29)  imply 

u£(x, -,p)  =  Ign(x, *,p)  +  0(p4)  as  p  0,  (4.31) 

uniformly  on  (-oo,tJ.  The  conclusion  under  case  (a)  is 
immediate  from  (4.31)  and  the  second  of  (4.30)  in  conjunction 
with  the  first  of  (4.1l). 

Next  consider  case  (b),  which  is  characterized  by  (4.7). 
Here  (4.6)  continues  to  hold,  whence  (4.31)  remains  valid.  In 
addition,  the  second  of  (4.7)  requires  the  three  first  moments 
(about  the  coordinate  planes)  of  the  tractions  on  each  family  of 
load  regions  to  vanish  up  to  time  t. 

Equations  (4.1), (4.2)  insure  that  fn  on  Dn  admits  the 
Taylor  expansion 

fn(a)  =  f n  +  in|pap  +  in(£)i  £  =  (a^ag),  (4.32)23 

where 

fn  h  fn(0)  =  ln  ,  =  1^(0)  (p=l,2).  (4.33) 

The  remainder  j>n£  C2(Dn)  and 

in(a)  =  0(a2)  as  a  s  |a|  0.  (4.34) 

Since  £  =  fn(a)  on  0n,  equations  (4.32), (4.7),  in  view  of  the 
first  of  (4.4),  at  once  furnish 

Recall  Footnotes  No.  19*22. 
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?J|p  |  S(i,T,p)ap(i)dA  =  -  J  Tl»J(a(l))S(l,*t,p)dA  (4.35) 

;n  JAn 

Ap  Ap 

for  all  ('tjp)  €  (-00  ,t]x(o,P1) .  On  estimating  the  right-hand 
member  of  (4.35)  with  the  aid  of  (4.34), (4.23),  and  (d)  in 
Definition  4.2,  one  infers  that 

fj|p  J  S(£,T,p)cp(£)dA  =  0(p4)  as  P  ->  0,  (4.36) 

A? 

uniformly  for  all  *t€(-oo,t],  Now,  (4.36)  may  be  regarded  as  a 
system  of  three  (inhomogeneous)  linear  algebraic  equations  in  the 
two  unknowns 


l  5 

in 


S(£,T,p)ap(£.)dA  (p=l,2) 


Furthermore,  because  of  the  last  of  (4.2),  the  coefficient -matrix 
of  this  system  has  the  rank  two.  Hence  (4.36)  imply 


S(l,-,p)ap(£)dA  =  0(p4)  as  p  — >  0  (p=l,2),  (4.37) 

In 
P 

uniformly  on  (-co,t).  From  (4.37)  and  the  second  of  (4.25),  in 
turn,  follows  the  estimate 

(4.38) 

uniformly  on  (-co,t].  Finally,  (4.38)  together  with  (4.31)  and 
the  first  of  (4.1l)  imply  the  first  of  (4.5)  with  6*4. 


l£n(x#*jP)  “  °(P4)  as  P  0, 
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Consider  at  last  case  (c).  Here  the  ordinary  equilibrium 

conditions  (4.8)  are  presumed  to  hold,  whence  the  loading  on  each 

family  of  load  regions  is  self -equilibrated .  In  addition,  as 

required  by  (4.9),  the  tractions  on  each  family  Ap  (n  fixed)  — 

up  to  time  t  —  are  now  supposed  to  form  a  parallel  system  at 

every  instant  and,  according  to  (4.10),  must  not  be  parallel  to 

oIl 

the  tangent  plane  of  B  at  £  . 

Since  (4.8)  include  the  assumption  (4.6)  underlying 
case  (a),  equation  (4.31)  is  satisfied  also  in  case  (c).  Further, 
equations  (4.8),  by  virtue  of  (4.4)  and  (4.32),  yield 


for  all  (t,p) € (-00 ,t]x(0,p1) .  Fran  (4.39)  one  draws25 

f^p  A  J  S(£,x,p)ap(£)dA  =  0(P4)  as  p  — >  0,  (4.40) 

uniformly  for  all  x€(-co,t).  Next,  substitute  from  (4.9)  into 
(4.40)  to  obtain 

ijpA^T)  j  p)ap(l)dA  ■  0(p4)  as  p  -»  0, 

Ap  (4.41) 

__ — - - 

It  follows  from  the  assumed  smoothness  of  the  boundary  B 

within  each  embedding  region  Hn  that  the  instantaneous  trac¬ 
tions  on  A^,  for  fixed  n  and  sufficiently  small  p,  cannot  at 
present  be  parallel  to  B  anywhere  within  Ap. 

25  Cf.  the  estimate  of  the  right-hand  member  of  (4.35). 


562(25)/21 


64 


uniformly  for  all  T€(-oo,t].  Scalar  multiplication  of  (4.4l) 

by  f1?  (q*l,2)  leads  to 

I  ° 

r  <Dn(l,T, p)Op(l)dA  =  0(p4)  as  p  -4  0  (p=l,2), 
Ap  (4.42) 

uniformly  for  all  T€(-co,t].  Now  note  from  (4.2)  that 
f^A  f^g  is  a  non-zero  vector  that  is  normal  to  the  boundary  B 

o  n 

at  v .  Consequently  and  by  the  hypotheses  on  k  ,  the  coefficient 

of  the  integral  in  (4.42)  is  uniformly  bounded  away  from  zero  for 

all  *r€(-ao,t],  whence  —  using  (4.9)  once  more  — 

f  S(£,-,p)a  (5)dA  *  0(p4)  asp-40  (p=l,2),  (4.43) 

Jn  p 

AP 

uniformly  on  (-ao,tj.  But  (4.43)  is  identical  with  (4.37)  and 
thus  the  first  of  (4.5)  with  6*4  follows  in  the  same  manner  as 
in  case  (b) . 

Wils  completes  the  proof  as  far  as  the  required  orders 
of  magnitude  of  the  displacements  are  concerned.  To  reach  the 
corresponding  conclusions  regarding  the  strains  one  may  proceed 
through  the  identical  argument,  taking  the  second  of  (4.11)  and 
the  second  integral  representation  in  (4.12)  as  the  point  of 
departure .  Since  the  foregoing  reasoning  in  no  way  depended  upon 
the  specific  nature  of  the  singularities  inherent  in  the  Green's 
displacement  01  entering  the  first  of  (4.12),  the  desired  con¬ 
clusions  for  the  strains  are  immediate  from  those  we  have 
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established  already.  Finally,  the  requisite  orders  of  magnitude 
of  the  stresses  follow  trivially  from  those  pertaining  to  the 
stratas  because  of  the  last  of  (4.1l)  and  the  third  of  (4.12). 

The  proof  of  Theorem  4.1  is  now  complete  in  its  entirety. 

Observe  that  the  content  of  Theorem  4.1  is  not  weakened 
if  the  conclusions  6-2,  6*3*  b=4  are  replaced  by  6  £  2,  6  2  3, 

6  >  4,  respectively.  'Hie  essential  significance  of  the  first 
conclusion  (6»2)  lies  in  the  fact  that  the  displacement^  strains, 
and  stresses  (at  a  fixed  interior  point  of  the  body)  are  bound 
to  vanish  at  least  to  the  order  0(p2)  as  p  ->•  0  in  the  absence 
of  any  restrictions  upon  the  loading  beyond  those  implied  by 
Definition  4.2.  In  particular,  thiB  order  of  magnitude  prevails 
regardless  of  whether  or  not  the  loading  on  each  family  of  load 
regions  is  or  is  not  self -equilibrated . 

In  case  (a),  where  (4.6)  is  met,  so  that  the  resultant 
force  belonging  to  each  family  of  load  regions  vanishes,  a 
reduction  of  the  (maximum)  order  of  magnitude  from  0(p2)  to  O(p^) 
is  guaranteed.  But  (4.6),  though  sufficient,  is  clearly  not 
necessary  for  6«3.  Analogous  comments  apply  to  the  further 

q  h 

reduction  from  0(P5)  to  0(p  )  assured  in  case  (b)  and  case  (c). 

Conditions  (4.7)  evidently  imply  the  ordinary  equilibrium 
conditions  (4.8).  The  converse  is  however  not  true,  whence  (4.7) 
represent  a  stronger  restriction  upon  the  loading  than  do  (4.8). 
On  the  basis  of  the  traditional  statement  of  Saint -Venant 's 
principle  in  elastostatics,  discussed  at  the  beginning  of  this 
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section,  one  would  expect  (4.8)  by  themselves  —  i.e.  In  the 
absence  of  the  additional  requirement  (4.9), (4.10)  that  the 
tractions  on  each  family  of  load  regions  be  parallel  and  not 
tangential  to  the  boundary  —  to  guarantee  the  reduction  to  6*4 
from  6=3  In  case  (a) .  That  this  expectation  Is  not  borne  out  by 
the  facts  Is  apparent  from  von  Mlses 1  [6]  counter-examples,  which 
refer  to  the  special  case  of  the  elastic  solid. 

Ihe  preceding  conclusion  has  a  counterpart  In  the  theory 
of  basic  singular  states  dealt  with  In  Section  2.  Thus  consider 
the  Kelvin-state  /°(x,t;0)  and  the  doublet-state  Ja^(x,t;0). 

Both  of  these  states  are  regular  for  all  (x,t) £ (E-R)x(-oo,oo ), 
where  R  is  an  arbitrary  (regular)  region  containing  the  origin, 
and  both  states  may  be  regarded  as  Induced  on  (E-R)x (-00,00 )  by 
their  respective  surface  tractions  on  the  boundary  II  of  E-R. 

Let  La,  Ma  and  LaP,  MaP  denote  the  resultant  force  and  the 
resultant  moment  about  the  origin  of  the  tractions  on  II  belonging 
to  ja  and  }a&,  respectively.  Then,  from  Theorem  2.2  and 
Theorem  2.3#  on  [0,oo), 


L°  /  0  ,  M®  =  0  , 

L®P  =  0  ,  /  0  (o^p), 

LaP  ■  0  ,  MaP  =  0  (a=P) . 


b  (*M4) 


Next  consider  the  rate  of  decay  of  the  corresponding  stresses  at 
infinity  (i.e.  "at  distances  large  compared  to  the  size  of  the 
region  of  load  application  II").  An  inspection  of  (2.27), (2. 35) 


562(25)/21 


67 


confirms  at  once  that  for  every  t  € [0,00 ),  as  x  s  |x|  — >  oo , 

Oa(x,t;0)  =  0(x"2),  Oa(x,t;0)  }  0(x"6)  (6  >  2), 

oaP(x,t;0)  =  0(x-3),  o°P(x,tjO)  /  0(x’6)  (6  >  3),  (a/p),M4.45) 

OaP(x,tjO)  =  0(x‘3),  oaP(x,t;0)  ^  0(x"6)  (6  >  3),  (a=p) . 


Consequently,  whereas  the  stresses  decay  more  rapidly  when  the 
resultant  force  of  the  loading  on  n  vanishes  than  when  this  is 
not  the  case,  the  additional  vanishing  of  the  resultant  moment 
fails  to  give  rise  to  a  further  reduction  in  the  order  of 
magnitude  of  the  stresses  as  x  — y  oo  . 

The  Saint-Venant  principle  contained  in  Theorem  4.1  may 
be  extended  to  accommodate  concentrated  surface  loads  with  the 
aid  of  a  corresponding  generalization2^  of  Theorem  3.3.  If  this 
is  done,  one  finds  that  Theorem  4.1  continues  to  hold  true 
provided  the  conclusions  6=2,  6=3,  6=4  are  replaced  by  6=0,  6=1, 
6=2,  respectively.  On  the  other  hand,  the  extension  of  the 
principle  to  anisotropic  viscoelastic  materials  would  require  an 
Integral  representation  for  the  solution  of  the  second  boundary- 
value  problem  appropriate  to  such  materials  and  analogous  to  that 
deduced  by  Fredholm  [17]  in  the  linear  equilibrium  theory  of 
anisotropic  elastic  solids. 

Acknowledgment .  The  authors  are  greatly  indebted  to  M.  E.  Gurtin, 
who  read  the  manuscript  and  supplied  numerous  most  helpful 
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- - 

Cf .  the  remark  on  concentrated  surface  loads  at  the  end  of 

Section  3. 
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